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UNIQUENESS AND PROPERTIES OF DISTRIBUTIONAL 
SOLUTIONS OF NONLOCAL EQUATIONS OF POROUS 

MEDIUM TYPE 

FELIX DEL TESO, J0RGEN ENDAL, AND ESPEN R. JAKOBSEN 


Abstract. We study the uniqueness, existence, and properties of bounded 
distributional solutions of the initial value problem for the anomalous diffu¬ 
sion equation dtu — = 0. Here can be any nonlocal symmetric 

degenerate elliptic operator including the fractional Laplacian and numerical 
discretizations of this operator. The function : R ^ R is only assumed 
to be continuous and nondecreasing. The class of equations include nonlocal 
(generalized) porous medium equations, fast diffusion equations, and Stefan 
problems. In addition to very general uniqueness and existence results, we 
obtain stability, L^-contraction, and a priori estimates. We also study local 
limits, continuous dependence, and properties and convergence of a numerical 
approximation of our equations. 


1. Introduction 


In this paper, we obtain uniqueness, existence, and various other properties 
for bounded distributional solutions of a class of possibly degenerate nonlinear 
anomalous diffusion equations of the form: 


(1.1) dtU - = 0 in Qt:=R^x(0,T) 

(1.2) u{x,0)=uo{x) on 


where u = u{x, t) is the solution and T > 0. The nonlinearity (p is an arbitrary con¬ 
tinuous nondecreasing function, while the anomalous or nonlocal diffusion operator 
is defined for any ip G as 

(1.3) Cf^[ip]{x)=f (ip{x + z) - ip{x) - z ■ Dip(x)li^i<i') dp(z), 


where D is the gradient, l| 2 |<i a characteristic function, and p a nonnegative 
symmertic possibly singular measure satisfying the Levy condition f |zpAl d/i(z) < 
oo. For the precise assumptions, we refer to Section 

The class of nonlocal diffusion operators we consider coincide with the gener¬ 
ators of the symmetric pure-jump Levy processes 0 [3 [39] like e.g. compound 
Poisson processes, CGMY processes in Finance, and symmetric s-stable processes. 
Included are the well-known fractional Laplacians —(—A)5 for s € (0,2) (where 
d/j,(z) = Cn^s for some Cff,s > 0 [M1I7])j along with degenerate operators, and 
surprisingly, numerical discretizations of these operators! 

In the language of 
On one hand 


, equation (1.1) is a generalized porous medium equation, 
since ip is only assumed to be continuous, the full range of porous 
medium and fast diffusion nonlinearities are included: p{r) = for m > 
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0. This is somehow optimal for power nonlinearities since if m < 0 (ultra fast 
diffusion), then not only uniqueness, but also existence may fail [12) . On the other 
hand, since (f is only assumed to be nondecreasing, it can be constant on sets of 


positive measure and then equation (1.1) is strongly degenerate. This case include 


Stefan type of problems, like e.g. when ci, C 2 , T > 0 and 


Lp{r) = 


C 2 r, r < 0, 

ci(r — T)+, r > 0. 


Many physical problems can be modelled by equations like (1.11. We mention flow 


in a porous medium of e.g. oil, gas, and groundwater, nonlinear heat transfer, and 
population dynamics. For more information and examples, we refer to Chapter 2 
and 21 in [IS] for local problems, and to gSlIMllTlIlSlIlZI for nonlocal problems. 

A key result in this paper is the uniqueness result for bounded distributional 
solutions of (1.11 and (1.2). Almost half of the paper is devoted to the proof of 
this result. Once we have it, we prove a general stability result, and then we 
obtain other properties like existence, L^-contraction, and many a priori estimates 
from more regular problems via approximation and compactness arguments. As 
straightforward applications of all of these estimates, we then obtain the following 
results: (i) Convergence as s —>■ 2“ of distributional solutions of 

(1.4) dtu + {-A)2ip{u) = 0 in Qt, 
to distributional solutions of the local equation 

(1.5) dfU — A(p(u) = 0 in Qt; 

(ii) continuous dependence in (to, s) S (0,oo) x (0,2] for the porous medium equa¬ 
tion of [37]) 

(1.6) -I-(—A)5ri|u|™“^ = 0 in Qt, 

including for the first time also the fast diffusion range; and (iii) convergence of 
semi-discrete numerical approximations of a class of equations including ([T3 (Cf. 
(|2.7|) and (|2.8|) in Section|2.2|). 


The uniqueness result is hard to prove because of our very general assumptions 
on the initial value problem combined with a very weak solution concept - merely 
bounded distributional solutions. This combination means that many classical tech¬ 
niques do not work: Fourier techniques are hard to apply because the problem is 
nonlinear and the Fourier symbol of could be merely a bounded function, en¬ 
ergy estimates do not imply uniqueness because (f is not strictly increasing, and 
L^-contraction arguments do not apply since we do not assume additional entropy 
conditions (cf. e.g. |S] for the local case), or equivalently, additional regularity 
in time as in m (see the uniqueness result for so-called strong solutions). The 
(weighted) L^-contraction argument for ordered solutions given in [T3] avoids these 
additional conditions, but it cannot be adapted here since it strongly depends on 
the equation being like (1.6) with 0 < to < 1 and s G (0,2). Finally, since our 


solutions are not assumed to have finite energy, the classical uniqueness argument 
of Oleinik [33] cannot be adapted either. We refer to (33] US] for the local case, 
and the uniqueness argument for so-called weak solutions in m for results in the 
nonlocal case. 


For the local equation (1.5), uniqueness for bounded distributional solution was 
proven by Brezis and Crandall in [18] under similar assumptions on ip and uq. 
Their argument is quite indirect and rely on a clever idea using resolvents and their 
integral representations (fundamental solutions). In this paper, we adapt such an 
approach to our nonlocal setting. But because of the generality of our diffusion 
operators, we cannot rely on explicit fundamental solutions for our proofs. Instead, 
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we have to develop this part of the theory from scratch, using the equation and the 
regularity that comes with our solutions concept to obtain the necessary estimates. 
To do this, a key tool is to approximate the possibly singular integral operator by 
a bounded integral operator and then carefully pass to the limit. This proceedure, 
and hence also the proof, is truly nonlocal - there is no similar approximation 
by local operators. The proof necessarily becomes much more involved than in 
[18j . and includes a number of approximations, a priori estimates, -contraction 
estimates, comparison principles, compactness and regularity arguments. It also 
includes new Stroock-Varoupolous inequalities and a new Liouville type of result 
for nonlocal operators. Both our approach and intermediate results should be of 
independent interest. 

Let us give the main references for the well-posedness of the Cauchy problems 
for (1.1) and (1.5). We start with the local case ( |1.5[ ). In the linear case, when 

When ip{u) = m™, it is a 
1. In the 


^p{u) = u, it is the classical heat equation, cf. e.g. 
porous medium equation, and a very complete theory can be found in 


general case, (1.5) is a generalized porous medium equation (or filtration equation). 
We refer again to [H]. Uniqueness of distributional solutions of this equation was 
proven in [18] for bounded initial data and continuous, nondecreasing (/?, and in [28] 
for locally integrable initial data, <f{r) = r™ for 0 < m < 1, and with regularity 
assumptions on dtu. Some nonuniqueness results can be found in e.g. [HI H5] . 
In the presence of convection, or if general L^-contraction results are sought, then 
so-called entropy solutions are a useful tool to obtain well-posedness mm- A 
very general well-posedness result which cover the case of merely continuous ip can 
then be found in [5]. 

In the nonlocal case, one linear special case of (1.1) is the fractional heat equa¬ 
tion dtu+ (—A) 5It = 0 for s G (0, 2). As in the local case, the initial value problem 
has a classical solution u{x,t) = *u{-,Q)){x) for It 

is well-posed even for measure data and solutions growing at infinity lini]. The 
fractional porous medium equations (1.6) are examples of nonlinear equations of 
the form (1.1). In 133 [37], existence, uniqueness and a priori estimates for ( |1.6| ) are 
proven for so-called weak L^-energy solutions - possibly unbounded solutions with 
finite energy. In m there are existence and uniqueness results for minimal distri¬ 
butional solutions of (1.6) with 0 < to < 1 in weighted L^-spaces (solutions can 
grow at infinity). We also mention that logarithmic diffusion {p>{u) = log(l -I-m)) is 
considered in [38] . singular or ultra fast diffusions in |12] . weighted equations with 
measure data in and problems on bounded domains in [ia[iniii7j. Energy 
solutions of equations with a larger class of nonlinearities ip and nonlocal opera¬ 
tors are studied in the recent paper [35] ■ The authors obtain results on well- 
posedness, continuity/regularity, and long time asymptotics. The setting, solution 
concept, and techniques are different from ours. Their operators can have some 
T-dependence, but the (singular part) must be comparable to a fractional Laplacian 
(i.e. be nondegenerate). Initial data in L°° n is assumed for uniquenss. In the 
T-independent case their assumptions are less general than ours, especially those 
for and the regularity of the solutions. Other types of equations of the form 
(l.I) can be found in jB). These equations involve bounded diffusion operators that 
can be represented by nonsingular integral operators of the form (1.3). Because of 
this, at least the well-posedness is easier to handle in this case. 

It should be clear from the previous discussion that even if our uniqueness result 
is very general, it is usually not strictly comparable to the other results. E.g. a price 
to pay to work with general ip and a very weak solution concept, is that solutions u 
have to be bounded. Our method of proof also requires that u — uq G L^{Qt)- For 
particular choices of ip, these assumptions may not be optimal. E.g. if you change 












F. DEL TESO, J. ENDAL, AND E. R. JAKOBSEN 


the solution concept and assume finite energy, then there are uniqueness results for 
unbounded solutions of ( |1.6[) in in |36l[37j. There are even uniquness results 
in weighted L^-spaces, see |l5j. Here the solutions are allowed to grow at infinity, 
but the uniqueness result is weaker in the sense that it only holds for minimal 
distributional solutions. 

There are other ways to generalize the porous medium equation to a nonlocal 
setting. In [n [m Sol [ini 141) the authors consider a so-called porous medium 
equations with fractional pressure. These equations are in a divergence form, and 
no uniqueness is known except when = 1. Finally, we mention that in the 
presence of (nonlinear) convection, additional entropy conditions are needed to 
have uniqueness as in the local case. Nonuniqueness of distributional solutions is 
proven in [2], and several well-posedness results for entropy solutions are given in 
[H [22l El] • These latter results requires ip to be linear or locally Lipschitz and hence 
do not apply to our case where Lp is merely continuous. 

Outline. In Section we state the assumptions and present and discuss our main 
results. The proof of the uniqueness result is given in Section]^ This proof requires 
a number of results and estimates for a resolvent equation - an auxiliary elliptic 
equation - and these are proven in Section]^ In Section]^ we prove the main stabil¬ 
ity and existence result, along with a number of a priori estimates. We then apply 
these results to prove the convergence to the local case, continuous dependence, 
and the properties and convergence of the numerical scheme in Section]^ Finally, 
after Section there is an appendix with the proofs of some technical results. 

Notation. For x S K, a;+ := maxjx, 0}, x~ := (—x)+, and sign''‘(x) is -1-1 for 
X > 0 and 0 for x < 0. We let B{x,r) = {y € : \x — y\ < r}, 1 ^( 2 ;) be 1 for 

X G H C and 0 otherwise, and supp ip be the support of a function ip. Derivatives 

are denoted by ', dt, and Dip and D'^ip denote the x-gradient and Hessian 
matrix of ip. Convolution is defined as / * y(x) = [f * g] (x) = /(x — y)g(y) dy, 

and (/, y) = fgdx whenever the integral is well-defined. If /,y G 
we write if,g)L^(M")- The L^-adjoint of an operator T is denoted by T*, and the 
reader may check that {£^)* = (see below for the definition of y*). A modulus 
of continuity is a nonnegative function A(e) which is continuous in e with A(0) = 0. 
By a classical solution, we mean a solution such that the equation holds pointwise 
everywhere. 

Function spaces: Cq, Cb, and Cp° are spaces of continuous functions that 
are vanishing at infinity; bounded; bounded with bounded derivatives of all orders; 
and smooth functions with compact support respectively. (^([O, T] ;Tioc(R^)) is the 
space of measurable functions ip : x [0,r] —>• K such that (i) ip{-,t) G 

for every t G [0, T]; (ii) for all compact K C \ip{x, t) — ip{x, s)| dx —t 0 when 

t-)’ s G [0,T]; and (in) ||^/’||c([o,T];Li(if)) := esssup(g[o_T] Jk dx < 00 . 

Measures: Sa(x) denotes the delta measure centered at a G K'^. Let X C 
be open and y a Borel measure on X. For x € X and H C A Borel, we denote 
ya;(H) = y(H -|- x) where il-|-x = {y-|-x:yG D}. Moreover, y* is defined as 
y*(H) = for all Borel sets B, and we say that y is symmetric if y* = y. 

The support of a Borel measure y on is 

supp y = {x G A : y(i?(x, r) n A) > 0 for all r > 0}. 

The Lebesgue measure of is denoted by drc if w is a generic variable on K^. 
Moreover, the tensor product dy(z) dw is a well-defined nonnegative Radon mea¬ 
sure since y is cr-finite (for more details, consult [3 Section 2.1.2].) 
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For the rest of the paper, we fix two families of mollifiers us, ps defined by 

(1.7) a;5(a) := 

for fixed 0 < w G satisfying suppw C i3(0,1), uj{a) = uj{—a), f uj = 1; 

and 

( 1 . 8 ) 

for fixed 0 < p G C'“([0,r]), suppp C [-1,1], p(t) = p(-t), f p = 1- 

2. The main results 

In this section, we present the main results: first of all uniqueness, and then 
stability, existence and a number of estimates for the solutions of O) and ( |1.2| ). As 
an application of our main results, we give compactness and continuous dependence 
estimates. We introduce a semi-discrete numerical scheme for even more general 
equations and show that convergence and other properties easily follow from our 
previous results. Finally, we establish a new existence result that also cover local 
diffusion equations. 

Throughout the paper we assume that 

(A^) 


(p : 


(A„p) uo G L° 


is continuous and nondecreasing; 


(A^) /i is a nonnegative symmetric Radon measure on \ {0} satisfying 


' d/r(z) 


1 dp{z) < oo. 


I\z\<l J\z\>l 

Remark 2.1. (a) Without loss of generality, we can assume (p(0) = 0 (by adding a 
constant to ip). 

(b) A nonlocal operator defined by (1.31 is a nonpositive operator (see Lemma 3.7). 


We use the following definition of distributional solutions of (1.11 and (1.2). 

Definition 2.2. Let uq G and u G L\aci.QT)- Then 

(a) u is a distributional solution of equation (1.1) if 


dtu - C^[p{u)] = 0 in V'{Qt), 


(b) u is a distributional solution of the initial condition (1.2) if 


ess lim 
(-> 0 + 


u{x,t)if{x,t)dx = / Uo{x)'tp{x,0) dx Wif G C° 


cA 


X [0,T)). 


The equation in part (a) is well-defined when e.g. (A,^) and (A^ I hold and u G 
L°°{Qt)- Note as well that the initial condition uq is assumed in the distributional 
sense {uq is a weak initial trace). See Lemma 2.21 below for an equivalent definition. 
We state the main result of this paper. 


Theorem 2.3. Assume (A,^) and (A,, |. Let u{x,t) andu{x,t) satisfy 
( 2 . 1 ) 

( 2 . 2 ) 

V'iQT) 


u,u G L°°{Qt), 
u-uG L^{Qt), 

(2.3) dtu - C^[p{u)] = dtu - C^[p{u)\ 

(2.4) ess lim /" {u{x,t) — u{x,t))fj{x,t) dx = Q for all tp G Cf° x [0,T)). 
Then u = u a. e. in Qt- 


in 
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Sections and 1^ are devoted to the (long) proof of this result. 


Corollary 2.4 (Uniqueness). Assume (A„q) and (A^). Then there is at 


most one distributional solution u of (1.1) and (1.2) such that u S L°°{Qt) and 
u- uo G L^{Qt)- 

Proof. Assume there are two solutions u and u. Then all assumptions of Theorem 
2.3 obviously hold (||u — u||li < ||u — mqIIli + 11^ ~ < oo), and u = u a.e. □ 


Remark 2.5. Uniqueness holds for uq ^ L^, for example uo(x) = c + 4>(x) for c € 


and (j) G n However, periodic uq are not included. In Section 2.3 

below we discuss some extensions of the uniqueness result. 

Next, we study under which assumptions solutions of 

(2.5) dtUn - C^’^[iPniun)] = 0 in Qt, 
converge to solutions of 

(2.6) dfU — C[ip{u)] = 0 in Qt- 



, Tn 
is a 


(i) -G Clf)] in Li(K^) for all ip G C)?°(K^); 

(ii) ipn ^ locally uniformly; 

(Hi) Un ^ u pointwise a.e. in Qt; 
then u is a distributional solution of ( |2.6[ ). 

This result is proven in Section]^ 


Remark 2.7. The limit operator C need not satisfy (A^ ), we can recover any oper¬ 
ator of the form C[ip] = D'^ip] + C^\ip]: the general form of the generator of a 

symmetric Levy process [7]. See sections 2.2 and 5.2 for more details and examples. 
An extension of this result will be discussed in Section [TS] below. 


The stability result will be used along with approximation and compactness 
arguments to obtain the following existence result and a priori estimates. 


Theorem 2.8 (Existence and uniqueness). Assume (A^), anduo G L°°(M^)n 

L\R^) 
isfying 


Then there exists a unique distributional solution u of (1.1) and ( |1.2[ ) sat- 
u G L°-{Qt) n L\Qt) nC{[0,T];Ll,{R^)). 


Remark 2.9. Existence results for merely bounded (and more general) initial data 
can be found in Theorem 3.1 in [15] in the setting of the fractional porous medium 
equation (1.6) with 0 < m < 1. 


Theorem 2.10 (A priori estimates). Assume { 


1, 

Am 

), uqjUq g l°°(r^) n 

L^{R^). Let u,u be the distributional solutions o 

f ( 

1.1 

with initial data uo,uo in 


the sense of Definition 2.2 (b), respectively. Then 

(a) (L^-contraction) (^(a:, t) —?).(a;, t))’*'da; < (uo(a;) —Uo(x))“''da;, t S [0, T]; 

(b) (Comparison principle) If uq < uq a.e. in R^, then u < u a.e. in Qt; 

(c) (L Abound) ||u(-,t)||ii(HJV) < ||uo||li(hjv), tG [0,r]; 

(d) (L°°-bound) ||u(-, t)|Uc=o(Riv) < ||uo|Uc=o(Riv), tG [0,r]; 
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(e) (Time regularity) For every t,s G [0,T] and compact set K C 

\\u{-,t) - u{-,s)\\l^k) < Ko (|i- s|®) + CK,ip,uo,ti (|i- + |i- s|) 

where Xuoi^) = inax|o -|<5 ||uo ~ '^o(' + o')IIli(r«)) I-^I the Lebesgue measure 
of K, and for some constant C independent of K, ip, uq, and p, 

= C'l^lf sup |(p(r)| + l) [ min{|2;|^, l}d^(z). 

^ kKIMollr.oo ^ J\z\>0 


(f) (Mass conservation) If, in addition, there exists L,S > 0 such that \p{r)\ < L|r| 
for |r| < 5, then 

/ u{x,t)dx= / Uo(a:)dx, t G [0,T]. 

7h« Jm" 

These results are proven in Section]^ 

Remark 2.11. The condition \(p{r)\ < L\r\ in Theorem 2.10 (f) is sharp in the 
following sense: If (p{r) = r™ for any to < 1, then there is = —(—A)® such 
that positive solutions u of (1.1) and (1.2) has extinction in finite time and hence 

f u ^ f Uq. Simply take N gN and s G (0,2) such that to < see [37] for 

the details. 

We now present several applications of the previous results. 

2.1. Application 1: Compactness, local limits, continuous dependence. 

We start by a compactness and convergence result for very general approximations 
of <ini) and ( |1.2| ). 

Theorem 2.12 (Compactness and convergence). Assume C : C'^{Qt) L^{Qt), 

/i„ satisfies ( |A^| ), pn and p satisfy (A^), and uo,n G L°°(M^) H for every 

n S N. Then if {un}n&f^ "is a sequence of distributional solutions of (2.5) with 
initial data {wo.nlnGN in the sense of Definition 2.2 (b), and 
(i) sup„ Jj^|>omin{|z|2,l}d^„(z) < oo; 

(a) sup„ ||Mo.n||L“(B«) < 

(Hi) ^ C[i}] in Li(R^) for all i) G C“(K^); 

(iv) pn ^ p locally uniformly; 

(v) Mo,n uo in Tio(,(M^). 

Then 

(a) there exist a subsequence {unjjjGN and auG C([0,T];Ljqj.(M^)) such that 
Un in C'([0,r];Lj(,^(M^)) as j ^ oo; 


(b) the limit u from part (a) is a distributional solution of (2.6) and (1.2). 


The proof can be found in Section |5.1| Using this result, we study the case 
= —(—A)5, s G (0,2). As expected, we find that solutions of the fractional 
equation (1.4) converge as s —>■ 2“ to the solution of the local equation (1.5). 
Then we obtain a new result about continuous dependence on (to, s) for the porous 
medium equation of m, that is, equation (|1.6|). 


Corollary 2.13. Assume (A,^) and uq G n L^(IR'^). 


(a) The distributional solution Us of (1.4) and ( |1.2| ), converges in C{[0,T] 

as s —>■ 2 to a function u, and u is a distributional solution of (1.5) and (|1.2|). 
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(b) Let Un and u be distributional solutions of (1.6) and (1.2) with {m,s) = 
[nirnSn) and {m,s) = {fh,s) respectively. If 


then Un 


(0, oo) X (0,2) 9 {mn,Sn) 
-^mn C{%T]-Ll^{R^)). 


(to, s ) G (0, oo) X (0, 2], 


The proof of this result can also be found in section 5.1 

Remark 2.14. When uq G L^(R^), the authors of [37] show continuous dependence 
in (^([O, T]; L^(M'^)) for ( |1.6[ ) and ( |1.2[ ) for (to, s) G ( i ^ When 

TO < , we are in the fast diffusion range and Corollary 

first continuous dependence result for this case. 


2.13 


(b) provides the 


2.2. Application 2: Numerical approximation, convergence, existence. 

Surprisingly, our class of operators is so wide that it contains a lot of its own 
numerical discretizations! It even contains common discretizations of local opera¬ 
tors as well. We illustrate this by giving one such discretization, a basic and very 
natural one, and then analyzing the resulting semidiscrete numerical method for 
(1.1), or rather (2.7). We prove that it satisfies many properties including conver¬ 
gence, and conclude a second and more general existence result. Consider 

(2.7) dtu - {L'^ + C^) [v3(u)] =0 in Qr, 

where is defined as before and is a possibly degenerate local operator 

;= tr[tTcr^Il^^(a;)] 


where cr = (tJi, ...., crp) G P G N, and (7^ G R^. Note that -|- is 

the generator of a symmetric Levy process, and conversely, any symmetric Levy 
processes has a generator like (cf. (3). Moreover, equation 0) and ( |1.5[ ) 

are special cases of (2.7) since cr and p, may be degenerate or even zero. 

For any h > 0, we approximate ( |2.7[ ) in the following way, 

( 2 . 8 ) 

where 

(2.9) 

( 2 . 10 ) 


dtUh - {Lf + Cf) [ip{uh)] =0 in Qt- 


^):=E 


i=l 


ij)(x + Gih) -\- ip{x — cfih) — 2ip[x) 


Cf[f)\{x) := ^ {ifix + Za) - f){x)) p{za + Rh) , 

a/O 

and Za = ha, a = (oi, ..., 0 ^^) G R/, = |[—1,1)'^. This is a finite difference 
approximation of L'^ and quadrature approximation of C^. 

Remark 2.15. (a) When a = Ci, a standard basis vector of then L®’ = and 
Lffifix) = ^ classical finite difference approximation. 


5.2 


and 


5.3 


_ 

(b) Both Lf and Lf are in form (1.3) and satisfy ( |A^[ ): cf. Lemma 

(c) L-V'(cr) = ELi <ylD^fj{x)a, = 

(d) C^^[ii]{x) = EaGZN E+r,. V’Ca; + ^) - V'(a;) d^(^) « ‘^fW\{x). 

(e) To avoid p{Rh) which may be infinite, we do not sum over a = 0 in Cf. 

We now show that the scheme has many good properties, including convergence. 

XP 

llliciuiuiip J-iajUyllL^^ ^ Ji-!. , 

, and h > 0. 


Proposition 2.16 (Properties of approximation). Assume (A,^), (A^ 
Mo, Wo e L°°{R^)nL^ 
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(a) (Existence and uniqueness) There exists a unique distributional solution Uh G 
L°-{QT)f^L\QT)r^C{[Q,T]-Ll^{W^)) of @ and (ig. 


(b) (LP-stable) ||Mii(-,t)||LP(B«) < l|wo|li~(Riv)||Mo|lii(HN), pG [l,oo], t G [0,T]. 

(c) (L^-consistent) For all tp G 

II M ~ MliLqRW) —>■ 0 as h^O^. 

(d) (Monotone) If uq < uq a.e. in , then Uh < Uh a.e. in Qt- 

(e) (Conservative) If in addition, there exists S,L > 0 such that |</3(r)| < L|r| for 
|r| < 5, then for all t G [0, T] 


/ Uh{x,t)dx= / uo{x)dx. 
Jr” 9r« 


Proposition 2.17 (Compactness of approximation). Assume (A^ 

n 


pNxP 


a G 


Uq G n L^(M^''), and h > 0. Then there is subsequence of dis¬ 

tributional solutions Ufi of ( |2.8[ ) and ( |1.2[ ) that converges in C([0, T]; ^[^^.(K^)) as 
/i — >■ O’*" to some function u. Moreover, u G ^““(Qt) C A^(Qt) CC'([0, T]; ^[^^.(M^)) 
and u is a distributional solution of and ( |1.2[ ). 

Note that Proposition 1 2.1 7| also provide a new existence result: 


Corollary 2.18 (Existence for (2.7|). Under the assumptions of Proposition 17/^ 
there exists a distributional solution u G L°°{Qt) H L^{Qt) C (^([O, T]; Lj'-^^(M '^ )) of 
(p^ and (iL^. 


In many cases we can combine the compactness result with uniqueness results 
for the limit equations, and hence obtain convergence for the approximation. 

Theorem 2.19 (Convergence of approximation). Under the assumptions of Propo¬ 
sition 2. I'll */ addition either a = 0 or fi = 0 and a = I (the iden¬ 


tity matrix), then the distributional solutions Uh of (2.8) and (1.2) converges in 
C([0,T];L(„,(K^)) as /i —>■ 0+ to the unique distributional solution u G L°°{Qt) C 
Li(QT)nC([0,T];Li„,(M^)) of (1^ and (p). 


The proofs will be given in Section |5.2[ 

Remark 2.20. (a) Our approximation is well-defined and converge for any prob¬ 


lem of the type (2.7), including strongly degenerate Stefan problems and fast 


(b) 


(c) 


diffusion equations. The scheme and convergence result thus cover cases that 
have not been considered before in the literature. For nonlocal problems of this 
type, there are very few results, and only for locally Lipschitz ip [151 p . 

To obtain a fully discrete numerical method, it remains to (i) restrict the 
method to some spacial grid and (ii) discretize also in time. Time discretiza¬ 


tion is easier and leads to a problem that no longer has the form (1.1); we will 


discuss it in a future work. Restriction to a spacial grid can always be done 


after a change of coordinate system: see Section 2.3 below. 

The existence result is a result where existence for problems involving nonlocal 
operators are exported to problems involving the “closure” of this class of 
operators - namely, operators of the form L'^ -\- . The proof is completely 

different from proofs based on nonlinear semigroup theory; see e.g. Chp. 10 in 

1551 . and 1571 . 


2.3. Remarks and extensions. 
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Alternative definition of distributional solutions. 

(1) A more compact form that we will use in the proofs is the following: 


Lemma 2.21. Assume (A^), ( Am„[), ( A^I and u € Then u is a 

distributional solution of ( 1 . 1 ) and ( 1 . 2 ) if and only if 

(u{x,t)dt'f{x,t) + tp{u{x,f))C^['ip{-,t)\{x'^AxAt+ / uo{x)'ijj{x,0) dx = 0 

V / .Ibn 


<0 -'B” 

for all if & x [0,r)). 


The easy and standard proof is omitted. 


About the initial conditions. 

(2) The solutions provided by Theorem |2.8| belong to T]; and hence 

satisfy the initial condition in the strong Lj^^-sense: For all compact K C , 



\u{x, t) — Uo(2:)| dx —>■ 0 


as t —>• 0. 


(3) If the initial conditions are satisfied in the strong Lj^^.-sense, then they are of 
course also satisfied in the distributional sense of Definition 12.21 


Extensions of the uniqueness result Corollary \2.4\ 

(4) With the same proof, we also get uniqueness for the initial value problem for 
the inhomogenenous equation 


dtu + C^[ip{u)] = g{x,f). 

(5) A close inspection of the proof reveals that we can replace continuity of ip in 


(A,p) by continuity at zero, Borel measurability, and (p{u) S L°°{Qt) (cf. [E]). 


Extensions of the stability result Theorem \2.()[ 

(6) When ipn is independent of n, we only need weak convergence of in (i): 

C^^[if] —> £[if] weakly in for all if € C^{Qt). 

Moreover, by considering subsequences we can replace (iii) by —>■ u in 

LJoc(Qt)- These observations follow by slight changes in the proof of Theo¬ 
rem [2]6] in Section m 

(7) A general condition for L^-weak convergence of [H]: There exist a G 
and a nonnegative Radon measure p, such that for all A G 

(i) sup„/|^|^omin{|zp,l}d/r„(z) < oo; 

(ii) /|z|<i dp„(z) -S' tr (crcr^A) -h /|^|<^ zAz'^ dp(z); 

(iii) /|z|>i <^Tn(z) dp(z). 

Here £ = ti[aa^D^] -|-£^: see [H] for a general discussion and more examples. 


Defining the scheme (2.8) on a grid. 


(8) By a coordinate transformation x = Ay^ L 
L^o _|_ £/i where Iq := 


can be transformed into 


■ / 

0 ■ 

0 

0 


cNxN 


. Up to permu- 
is orthonormal, 


/ is an identity matrix, and dfl{z) = dp{A~^z) satisfies (A^ 
tations of the components oi y, A = QJ where Q G 
Qaa^Q^ = diag(Ai) for Ai > 0, and J = diag(^/ci) where Ci = 1 if Ai = 0 and 
Ci = ^ if Ai > 0 for z = 1,..., iV. 
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(9) For the new operator our approximations produce an operator 

that can be restricted to the (j/-)grid Qu := {h > 0 ), that is : 

—>■ is well-defined. 

3. The proof of uniqueness 

3.1. Preliminary results. A crucial part in the proof is played by the following 
linear elliptic equation 

(3.1) eve{x) - C^[ve]{x) = g{x) in 


where e > 0 and defined by (1.3). Its solutions will be denoted by 


B^[9\{x) ■■= Ve{x). 

Formally, = {el — is the resolvent of . Note that may be very 

degenerate and therefore Fourier techniques do not easily apply (cf. Example 3.1 
and Remark 3.8 (a) below). The main results about equation are given below, 
while most of the proofs will be given in Section Note that in [T5] such results 
are easy in view of an explicit representation formula for B^. Here, on the other 
hand, they are not easy and we have to work quite a lot to prove these estimates. 
The method of proof is different, more nonlocal, and requires less of the operator. 


Theorem 3.1 (Classical and distributional solutions). Assume (A^) and e > 0. 

(a) If g € C“(K^), then there exists a unique classical solution B^[g] G 
of (3.1). Moreover, for each multiindex a G 


(b) If g G then there exists a unique distributional solution B^[g] 

L\R^) of ( 1 ^. Moreover, 

(e) If g G then there exists a unique distributional solution B^[g] 

L-(M^) of ([rg. Moreover, 


emmi 


< 


Remark 3.2. If g e C L°°, then £||i3^[(7]||LP < HsHlSo || 5 ||£i for any p G (1, oo). 
When a smooth g depends also on time, then B^[g] will be smooth in time and 


space. 


Corollary 3.3. Assume (A^), e > 0, and 7 S C, 


X [0,T)). Then 


(a) B^[j] G C£°(M^ X [0,T)). 

(b) B^['y\{x,-) is compactly supported in [0,r). 

(e) dt {B^ [ 7 ]) = B^ [da] and B^ [ 7 ], B^ [da], [B^ ^GL^Qt). 

Proof, (a) A standard argument using difference quotients, linearity and uniqueness 
of the problem, the L°°-bound of Theorem 3.1 (a), and induction on n, gives that 

(3.2) a£L'“H£[ 7 ] = B^[d{fD°‘j] in Qt 

for every n S N and a G N^. This argument is almost exactly the same as the one 


given in the proof of Proposition 6.8 (d) below. Then by Theorem 3.1 (a), 

el]dpn^B,^[j]h^(Qa < ]]drB^jh^(Qa. 

(b) Holds since B^ is an operator in the spatial variable x and B^[0] = 0. 


(c) Note that dtB^^y] = B^[da] by (3.2), and by Theorem |3.1| (b) and the time 
continuity of 7 and B^[y], 

4B^M]LHQa < II7||li(Qt): 
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which is finite because 7 S C^{Qt)- Hence it follows that 

4dtiB^W\\LHQT) = 4B49tl]\\L^QT) < Wdah^QT), 

By equation ( |3.1[ ), = eB^lj] — 7 for all {x,t) S Qt- Since both - 8 ^( 7 ] 

and 7 are in ITJQx), it follows that also £'^[8^[ 7 ]] € L^{Qt)- 

The operator 8^ is self-adjoint in the following sense: 

Lemma 3.4. Assume (A^ ), g G f G and £ > 0. Then 


□ 


/ B^[g]{x)f{x)dx = / g{x)B^[f]{x)dx. 

Jr" Jr" 

The proof is given in section To prove these and other results in this paper, 
we will need some properties of the nonlocal operator that are given below. 


Lemma 3.5. Assume (A^ 
(a) If^GC^ 


then 


\C^^[4\{x)\ < - max \D'^^{x + z)\ 
Z |z|<l 


ld<i 


\z\^d^l{z) 


'ld>i 


dn{z). 


(h) Let p G { 1 , 00 } be fixed. If ip G then 


\m4\\ 


LP(R") ^ 


< 4\D^ 


|LP(R«) 


/ NPdAi(z)-b 2||'!/'||LP(Riv) / dp{z). 
J\z\<i •'bl>i 


(c) If ipi G VL2,i(KAf) andiP 2 e then 

[ ilJiC^[ip2]dx = [ C^[ifji]ijj2dx. 

Jr" Jr" 


Remark 3.6. (a) If '0 G n then C^[ijj]{x) is well-defined by (a), 

(b) If < 00 , a density argument and the symmetry of p reveals that 


C^[4>]{x) = [ U{x + z)- cpix)) dp{z), 
J\z\>0 ^ ^ 


and the assumptions of Lemma 3.4 can be relaxed to g G ), / G 

for p G { 1 , 00 }, and 0i G L^{W'‘) and 02 G respectively in (a), (b), 

and (c). The second derivative part of the estimates in (a) and (b) then have 
to be dropped and the remaining term modified accordingly. 


A proof of Lemma 3.5 can be found e.g. in Sections 1 and 4 in [3]. 

and Ip G Cpfi{R^). Then 


Lemma 3.7. Assume (A, 




where 

(xcp(i) 

Moreover, crcp-iO > 0 and 





L2(RN) 


cos(z • fP) dp{z). 




2 

L2(KN) ■ 


Remark 3.8. (a) (T£p is the Fourier symbol of . In our generality it may not be 
invertible or have any smoothing properties. An extreme example is p, = Szg for 
Zq fiz 0, where ocpifif) = 1 — cos^o ’ 4 this is a bounded function with infinitly 
many zeros. 

(b) If 0,£^[0] G then a density argument shows that the Fourier symbol 

exists and the conclusions of Lemma [3.71 still hold. 
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(c) The notation 2 is used to denote the square root of the operator in the 
Fourier transform sense. 

Proof. By the definition of Fubini’s theorem, and the symmetry of /r, 

? f 




/r« 


'| z |>0 


‘ + z) — — z ■ D'i/)(a;)l| 2 |<i^ d^(z) dx 

• Cl|z|<i-^(V')(0) dn{z) 
= -^(V')(0 / (cos(0-C) - 1 ) d/i(z). 


'| 2|>0 


'| 2:|>0 

To show t he se cond part of the lemma, note that > 0 and if, G 
(cf. Lemma 3.5 (b)). It follows that P{ijj), aci^P{ip) € and then by the 

inequality 2ab < + 6 ^, T{ijj) € By Plancherel’s theorem, 

= = (rw, -„,,r 

= - = - 

which completes the proof. 

The following theorem is a key technical tool in our uniqueness argument. 




L2(KJ' 


[L^Hif] 


L 2 (R«) 


□ 


Theorem 3.9. Assume (A^l and supp/x ^0. If v G C'o(ffi^) solves 

=0 in V'{R^), 

then V = 0 for all x € K^. 

We give the proof of Theorem 3.9 in Appendix]^ In the local case [TB] such a 
result follows for example from the Liouville theorem for the Laplacian. On one 
hand, our result is much weaker since we need to ask for some kind of decay at 
infinity. On the other hand, Theorem |3.9| covers very degenerate operators 
which do not satisfy any sort of Liouville theorem. 

Example 3.1. Let /i = S 2 Tr + f- 2 Tr- Note that ( |A^[ ) holds and that for smooth 
functions v, 

= v(x + 27r) — 2v(x) + v(x — 27r). 

The function v = cos G is an example of a nonconstant function that 

satisfies C^^\v\{x) = 0 in K, and hence the Liouville theorem does not hold for . 

3.2. The proof of Theorem |2.3[ We define 

U{x,t) := u{x,t) — u(x,t) and ^{x,t) := (p(u{x,t)) — ip{u{x,t)). 


By the assumptions ( |2.1[ ), ( |2.2| ), and (A,^, 

U G L^{Qt) n L°°{Qt), ‘f’ G L°°{Qt), 


and by (|2.3|), (|2.4|), and Lemma 2.21 
(3.3) 


/ 


Udtif + ^C^[if ]) da; dt = 0 for all if G C[ 


'oo 


X [0,T)). 


We emphasize that this equation also incorporates a zero intitial condition for U. 
We now define the function hg{t) which will play the main role in the proof: 

(3.4) hS)-={B^m,t),U{;t))=[ B^[U{;t)]{x)U{x,t)dx. 
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Note that G L^{{0,T)) since \\hs\\iM(o,T}} < 1 \\U\\l’^(Qj.)\\U\\l^Qt} by Theorem 


3.1 (b). For the proof of Theorem |2.3[ we will now show that there is a sequence 
0 + such that lime^_^o+ = 0 - To do that we start by the following 


lemma: 


Lemma 3.10. Assume (A^), U G L^{Qt) L°°{Qt), $ G L°°{Qt), and (3.3) 
holds. Then 

(a) jj (BnN]atV' + (e-Be^[^]-^)V')dxdt = 0 for all ip e x [0,T)). 

(b) Bf[U{-A)]ix) = J s)](x) — $(a;, s)^ ds a.e. (a;, t) G x ( 0 , T). 

(c) For a.e. t G (0,T), t)||ioo(RiV) < 2t||d>||icx=(Q^). 


Proof, (a) We fix 7 G C; 


'oo (mN 


X [0,T)) and take tp = as a test function 


in (3.3). Note that ip is an admissible test function by a density argument using 


Corollary 3.3 (a)-(c) and [/, <& G L°°{Qt). Then by ( |3.1| ) and Corollary 3.3 (c), 

0 = 11^ {udtm-f]) + [B^M ) dxdt 

= JJq (yB^[da] + d>(sB^[-/]-^'^^dxdt. 

Finally, the self-adjointness of B^ (cf. Lemma 3.4) yields 

^ (^B^[U]da + (eB^m - $)7) dxdt = 0, 

which completes the proof. 

(b) This result follows from (a) and a special choice of test function. For 0 < s < T, 
a > 0 , and 0 < i5 < T — a, we define 


1 


t < s — a 


Oait) = { I — — s + a) s — a < t < s 

0 t>s 


and 0a,sit) = 0a* ps{t), 


where the mollifier ps is defined in ( |1.8| ). Then 9a^s G (^“((OjT)) C L^((0,T)) 
and suppjda^i} C [—oo,T). Let 7 G and take ipix^t) = 9a^s{t)"f{x) G 

X [0, T)) as a test function in part (a). Then we use properties of mollihers 
and Lebesgue’s dominated convergence theorem to send i5 —>■ O'*' and get 

11^ [b^[UK + - d>)9a)j dxdt = 0. 

By Fubini’s theorem and since 9'^{t) = —^'ds-a<t<s and suppjda} = [0, s], we hnd 
that ^ ^ 

f f B^[U]dt+ f (eB^^[4>]-$)6»,,dt )7 dx = 0. 

7r« \a Js-a Jo J 

We now send a 0+. Since B^[U{-,t)]{x)"f{x) dx G L^(0,T) by Fubini’s 
theorem, 

- f f B^[U{-,t)]{x)^{x)dxdt ^ f B^[U{-,s)]{x)j{x)dx as a —>■ 0''' 

a Js-a Jr^ Jr^ 

for a.e. s by Lebesgue’s differentiation theorem. For the other term, we may use 
Lebesgue’s dominated convergence theorem to pass to the limit. Since 9a —> l[o.s) 
pointwise, we hnd that for a.e. s G [0,T], 

B^[U{-,s)]{x) + J (^eB^[4>(-,t)](x) - dt^ 7 (x) dx = 0. 
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Since 7 e is arbitrary, part (b) follows, 

(c) By part (b) and Theorem 


3.1 


(c), ||S|‘[17](-,t)||i=o(Riv) < 2t||$|ii=o(Q,^) a.e. □ 


Proposition 3.11. Assume ( |A^[ )^ U G L^{Qt) H ^ G L^{Qt), and 

(3.3) holds. Then h^{t) defined by (3.4) is absolutely continuous and 

K{t) = 2{sB^m;t)-^;t),U{;t)) in l?'((0,r)). 

The proof below is an adaptation of the proof in Ha pp. 157-158]. 


Proof. Let the mollifier ps = ps{t) be defined in (1.8), the extension U he U on Qt 
and zero outside Qt, and 

Us{x,t) := U{x,-) * ps{t) = / U{x,s)ps{t- s)ds. 

Jr 

By Young’s inequality, \\Us\\L--iQT) < I|C^IIl“(^ and \\Us\\l^Qt) < \\U\\lPQt)- 
Moreover, the time continuity of Us, Corollary |3.3| (c), and Lemma 3.4 yields 

d 


(3.5) 


dt , 




B^[Us]Usdx = 2 


dt {B^ [Us] )Usdx = 2 dt{Us)B^[Us] dx 


for t e K. 

Let us show that 


B^[Us{-,t)]{x) = / B^[U{-,s)]{x)ps{t - s)ds in Qt- 


(3.6) 


First assume that U S Cff(QT) H L^{Qt). Then B^\U[-,t)\ G C'((°(K'^) C 
for t G [0,T], and thus, it solves ( |3.1| ) pointwise in Multiply this equation 
by ps{s — t), integrate over K, and use Fubini’s theorem and the uniqueness in 
Theorem 3.1 (b) and (c) to find that (3.6) holds. A density/mollific ation argument 
using uniqueness and L^(IR^) and ) estimates from Theorem 3.1 then shows 

that (3.6) also holds (a.e.!) for U G L^(Qt) H L°°{Qt). _ 


Let the extension $ be $ on Qt and zero outside Qt- Using Lemma 3.10 (a) 
with test functions if G C'()°(]R^ x {S,T — <5)) we get that 


dtB^[Usi-,t)]ix)= i^{eB^[^-mx,-)*psj{t) a.e. in R^x{S,T-S). 

For any 0 G C'(?°((0,T)) and sufficiently small S, we then conclude from (|3.5|) that 


{B^[Us]{;t),Us{;t))Q'{t)dt = 2 r {{sB^m - ^)*ps{t),Us{-,t))e{s)dt. 

3 ^0 


Qt, 


By properties of mollifiers and Theorem |3.1| (b) and (c), 

Us ^ U in L^{Qt), 

{eB^[^] — ^) * ps ^ cBg [4>] — < 1 ) a.e. in 

e\\Bf‘[Us]\\L°°{QT) < II^IIl“(Qt); 

|(eBn$]-l>)*P5| <2||$||i»(Q^). 

Now we send (5 —)■ 0+ using Lebesgue’s dominated convergence theorem, and then 
by the definition of hg, we find that 

- r hS)Q'{t)dt = 2 r {eB^^[m-U)-H;t).U{;t))Q[t)dt. 

Jo Jo 

That is, he is weakly differentiable and the weak derivative is 
K{t) = - H;t),u{;t)). 
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Moreover, h'^ S L^{{0,T)) since by Theorem 3.1 (c), 

r \K{t)\dt<m\\L^iQ^)\\UhHQ^)■ 

Jo 

Hence, hi;{t) is absolutely continuous, and the proof is complete. 


□ 


Proposition 3.12. Assume Q, U G L\Qt) n L°°{Qt), $ G L°°{Qt) 

and (3.3) holds. Then 
(a) For a.e. t G [0, T] 

hs) = emm-Mh + 

(h) If a sequence —>■ 0 a.e. in Qt as e„ — >■ O'*', then for a.e. t G [0,r], 

lim heAt)=G. 

En—>-0+ 

We need a technical lemma (cf. [E]). 


Lemma 3.13. Assume ( |Ay[ ) and (2.2). Then the Lebesgue measure of the set 
:= {{x,t) G Qt ■ \(p{u{x,t)) - q}{u{x,t))\ > ^}, 
is finite for all ^ > 0. 


Proof. Define the set 

= ii^A) € Qt ■ \u{x,t) - u{x,t)\ > d}. 

If {x,t) G S^, then by the continuity of (p there exists a d > 0 such that \u{x,t) — 
u{x,t)\ > 6, that is, C S^. By (2.2), 

< // \u{x,t) — u{x,t)\dxdt < oo, 

JJ Qt 

and thus, also has finite Lebesgue measure. □ 


Proof of Proposition S.lt^ (a) By the assumptions. Theorem 3.1 (b) and (c), in¬ 
terpolation between and and Fubini’s theorem, we have for a.e. 

t G [0,T] that U,B^[U] G and 


(3.7) eB^[U]- Cf^[B^[U]] = U in V'{R^). 

Hence it follows that C^[B^[U]] G where is defined through the relation 

[ C^[B^[U]]ifdxdt= [ B^[U]C^^[fj]dxdt for all fjeC^{R^). 

Jm"' Jm"' 

Using Plancherel’s theorem and Lemma |3.7[ we then find that for any G 

[ F{C^[B^[U]])F{i:)df= [ FiB^[U])FiC>^[i;])df 

Jr'^ ./r« 


F{B^[U])ac.if)AWdf, 


and hence 

[ Fim)(B{c^mu]]){f)+ac.{OBmu]m) df = o. 

Jr" ^ ' 

Then by a density argument, we conclude that 

Ti^\B>f\U\\){i) = -ac.{f)T{B>f\U\m in L\R^), 

and thus, for a.e. t G [0,T], we have C^[B^[U]] = in L'^{R^). 
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Since U,B^[U],C^^[B^[U]] S equation ^ holds in L^(R^). By 

Lemma [3.71 Remark |3.8| (b), and the definition of (see (3.4)), we have for a.e. 
t G [0, T] that 

= (B^m-J),eB^lU](-,t) - 

=e\mu]i;t)\\hiMn + wicnHBnumh^^.y 


(b) By part (a), Proposition 3.11 and = {u — u){(f{u) — (f{u)) > 0, 

rt 


(3.8) 


0 < hs{t) = hs{0+) + f h'g{s)ds 

Jo 

< h,{0+) + 2 [ {eB^m-, s),U{; s)) ds. 
Jo 


By the (absolute) continuity of he, Holder’s inequality, Lemma 3.10 (c), and 
Lebesgue’s dominated convergence theorem (valid since U G L^{Qt)), 

he(0+)= lim y f heis)ds< lim ^ f ||R^[t/](-, s)||i=o(Riv)||17(-, s)||ii(Rjv) ds 
t^o+ t Jo ‘^0+ t Jo 

< 2||$||L»(Q,r) 1™ / l|C^(-,s)llLi(B«)l( 0 .t)(s)ds = 0. 

t^0+ Jq 


Let ^ > 0. By the self-adjointness of B^ (cf. Lemma [3.4[ ) and Theorem |3.1| (b), 
we get for a.e. t G [0, T] 


{eB^m;t),U{;t)) = [ d>{x,t)eB^[U{;t)]{x)dx 


< ll^ll 


< ll^ll 




\eB^[U] \ dT + C 


/{|$(a:,t)|< 5 } 


\eB^[U]\ dx 


/B« 


\eB^[U{-,t)]\ l|3.(^,t)|>4 dx-hC||H(-,<)||ii(RiV). 


Let t be a point where this inequality holds and enB^[U{-,t)] -G 0 a.e. x and 
|ei3^[?7(-, t)](a;)| < \\U\\l,x,(^q^) a.e. x (using Theorem 3.1 (c)). For any r? > 0, take 
^ such that ^\\U{-,t)\\i^i < J-rj. Then note that |£:i3^[{7]| l|$(a;_t)|>j is dominated 
by ||C^||L“l|$(x,t)|>{ which is integrable by Lemma 3.13 By Lebesgue’s dominated 


convergence theorem it then follows that |£„i?^^[{7(-,t)]|l|$(a;^t)|>j dx < 
when e„ is small enough. Since this holds for a.e. t G [0,T], we have proven that 

lim (£„R^^[<l)](-,t),{7(-,t)) < 0 for a.e. te[0,T]. 

>-0 + 

We conclude the proof using Lebesgue’s dominated convergence theorem to send 
£n —>■ 0+ in ( |3.8| ) (the integrand is dominated by ||<i)||ic=o(Q,j,)||{7(-,t)||ii(Riv) G 
L^{{0,T)) since U G (Qt) n L°°(Qt))■ □ 


Proposition 3.14. Assume ( |A^[ ), supp/x ^ 0, and g G H Then 

there exists a sequence such that enB^^[g] —)■ 0 a.e. in M.^ as e„ — )■ 0+. 

This proposition will be proven later in this section. We are now ready to prove 
our main result. 

Proof of Theorem\2.S\ In the case that supp^ = 0, ^ = 0 and = 0. Then equa¬ 


tion (1.1) becomes the ODE ut = 0, and uniqueness follows by standard arguments 
(e.g. one can easily deduce that |zi(x, t)— m(x, <)| dx < 111(3^) 0) — u(x, 0)| dx). 
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Now consider the case supp^ ^ 0. By Proposition 3.14 and 3.12 (a) and (b), 
there is a sequence such that for a.e. t G [0,T], 

(3.9) + ||(/:^)5[<[t/]](.,t)||i. ^0 as ^ 0+. 

Let i/^ € By Plancherel’s theorem, Lemma 3.7 and Cauchy-Schwarz’ 


inequality, and finally, by (3.9), we get for a.e. t € [0,T] that 










as e„ —>■ O'*'. Moreover, by Cauchy-Schwarz’ inequality and (3.9), we have for a.e. 

tG [0,T] 


'R» 


enB^^\U]'ip dx 


< ||e„B,^J[/]|U2(R»; 


L2(RiV) 


0 as e„ 


0 +. 


Hence we conclude that as —>■ 0+, 

U = enB^^[U]-C^[B^JU]]^ 
for a.e. t G [0,r]. That is, 

u — u = U = 0 in 


in 


V'( 


oAf 


for a.e. t G [0,r], and then a.e. in Qt by du Bois-Reymond’s lemma. 


□ 


In the rest of this section, we prove Proposition 3.14 For 7 G we 

let Vg := sB^[y] be the unique smooth classical solution (see Theorem |3.1| (a) and 
Corollary 3.3 (a)) of 

(3.10) eve{x) — C^[ve]{x) = £ 7 (x) for all x G . 

We want to prove that there exists a sequence such that = e„H^^[ 7 ] —0 as 
En —>■ 0+ for every x G and every 7 G C^{R^) . 

l^). Then there exists a sequence 


Lemma 3.15. Assume (A^ 


and 7 G (7° 


that converges locally uniformly in k" as £n 


]>N ^ ^ O’*'. Moreover, the 

corresponding limit v is uniformly continuous, lim| 3 ,|_^Qo 1 ; = 0 and satisfies 

/:^H(a:) = 0 m V'{R^). 

Lemma 3.16 (Barbalat). If if G L^{R^) is uniformly continuous, then 
lim if{x) = 0 . 

|a;|—>-oo 

For a proof, see e.g. Lemma 5.2 in |30] (take G = R^ and B = M). 


Proof of Lemma \3.15\ We recall that Ve := ei?^[ 7 ]. By Theorem 3.1 (a), 

for each multiindex a G N'^. So, then any sequence is equibounded and 

equilipschitz. By Arzela-Ascoli’s theorem, there exists a subsequence such that 
-G V locally uniformly as n —)■ cx). Since is uniformly continuous (the 
derivative of exists and is bounded) and by the local uniform convergence, for 
every ij > 0 and i? > 0 we can find some n > 0 such that max{|r;(a;) —W£^(a;)| : |a::| < 
R} < r]. Thus, we have the following estimate for every R> 0 and |a;|, \y\ < R, 

\v{x) - v{y)\ < |w(a;) - We„(a;)| -b \ve^{x) - Ve„{y)\ + \ve„{y) - v{y)\ 

<2r?+ ||L>7||ioo(Riv)|x-y| 

As R is arbitrary, v is Lipschitz continuous with Lipschitz constant ||L17 || l°°(r -”)i 
and thus, uniformly continuous. Furthermore, Fatou’s lemma and Theorem 3.1 (b) 
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give that ||T||ii < liminf„_j.oo \\veJ\L^ < ||7lUi- % Lemma 

0 . 


3.16 


lim 


|a;|—>-oo 


j{x) = 


Multiplying (3.10 1 by a test function, integrating over and using self-adjointness 
(cf. Lemma 3.51 of we get 

En [ ( Ve^C^[^p]dx = En [ 'ytjj dx for all tp G 

is™ Jm" Jr" 


Since llue 


< 


licxj by Theorem 3.1 (c), we use Lebesgue’s dominated conver¬ 


gence theorem to take the limit as —?> O’*" , to find that 

0 = lim f Vs^C^[ip]dx = f vC^[^p]dx for all ipGCp°(R^), 

e.x-lO+ Jrn Jrn 


which completes the proof. 


□ 


Lemma 3.17. Assume (A^ 
sequence [( 7 ]}„gN that converges in 


and g G 


n L°° 

as £. 


Then there exists a 
0 +. 


Proof. Note that Ug := £Bp-[g] is the unique distributional solution (see Theorem 


i>N\ 


3.1 (b) and (c)) of the following elliptic problem 

£Ue{x) — C^[Ue]{x) = £g{x) in !)'( 

(b) and (c) and the linearity of the above equation, for any h G , 
D Ikellii < IIsIIli and \\ue{- + h) - u^Wli < \\g{- + h) - g\\Li . 


3.1 


By Theorem 

Wueh^ < Ml 


3.1 


Now let K C be any compact set, and define wf[x) = Ue{x)lK{x). The 
uniform in £ bound ensures that the family M := {w^ }£>o C is uniformly 

bounded in Moreover, by continuity of the L^-translation, Theorem 

(b) and (c), and Lebesgue’s dominated convergence theorem, 

\\w^{- + h)-w^\\Li 

< II (We(’ + h) - Ue) 1 k{- + ^)||li + ll^e (lif (' + h) - Ik) IIli 

< II5(- + ^)-5 ||li+ II5||l” / \1k{x + h) - 1k{x)\ dx ^ d as |/i|0. 

JR" 

Combining the above results, we see that M is relatively compact by Kolmogorov’s 
compactness theorem (see e.g. [551 Theorem A.5]). Hence, there is a convergent 
subsequence in L^{K). 

Now, cover by a countable number of balls i?„. Then the above argument 
holds for K := Bn for every n € N. A diagonal argument then allows us to pick a 
subsequence which converges in L^{Bn) for each n, and thus in D 


Remark 3.18. By Theorem 3.1 (a) and Arzela-Ascoli, we can have D°‘Vs —>■ Wa 
locally uniformly in as £ —)■ 0+ for all multiindex a G N^. However, because 
of the lack of uniqueness in C^[v\{x) = 0, we do not know if D‘^v = Wa- Hence, we 
are forced to work with distributional solutions of C^[v]{x) = 0. 


Lemma 3.19. Assume (A^), g G L^{R^)nL°°(R^), and {£nBM9]}neN converges 
in V ^nBMl]{x) —0 as £„ —?► 0+ for every x G and every 7 G 

then EnB^g] —>■ 0 m as £„ —>■ 0 +. 


Proof. By the self-adjointness given in Lemma 3.4 and the definitions := 
£nBp[g], Vsn '■= SuBplj], we have 
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Since WvsJIl^^ < UtIU” by Theorem [3l] (c), \g{x)ve^{x)\ < \g{x 
the assumption and Lebesgue’s dominated convergence theorem, 


i=o. Then by 


lim f Ue^{x)^{x) dx = 0 for all 7 G 
£^->■ 0 + J-^N 

Hence > 0 in and since the distributional and limits coincide (by 

uniqueness), it follows that Ue„ —t 0 in as £„ —>■ 0 +. □ 


Proof of Proposition \3 . lf\ Let 7 G be arbitrary, and recall the definitions 

= Ve and eB^[g] = Ug. Lemma [3.15| yields a subsequence such that —>■ v 

locally uniformly as e„ —>■ 0+ with v G C'o(K^) and C^[v]{x) = 0 in _D'(K^). Then, 
Theorem 3.9 ensures that v{x) = 0 for every x G 


Hence, Lemma 


3.17 


and 


3.19 


give that —>• 0 in as £„ —>■ 0+. Finally, 

take a further subsequence (still denoted by £„) such that —>■ 0 a.e. in as 

£„ 0 +. □ 


4. Stability, existence and a priori results 

In this section, we will start by showing the stability result stated in Section 
and then we continue by showing existence and a priori results for ( |1.1[ ) . The latter 
part will follow by regularization and compactness from results in |23| for the case 
ip G and Mo e L°°(K^) n 


Proof of Theorem 2.6 Since Un are distributional solutions of (1.1), we will take 
the limit as n —)■ 00 to see that so are also u. 

Assumption (iii) and the uniformly boundedness of ||m„||loo(q,j,) gives for all 
Ip G C^{Qt) that 


10 JR'^ 


Undtip dx dt 


10 JR" 


udt^pdxdt as n —)■ 00 . 


To prove convergence of the £^’*-term in the distributional formulation we pro¬ 
ceed as follows 


/o Jr" 

rT 


(pPn{un)jO.^"[pj] - dxdt 


/ / ipniUn){C^"[tp]- C[lp])dxdt+ / {tpniun) - ifiun)) C[lp] dx dt 

IQ JR" Jo JR" 


10 JR" 


{(p{Un) - ip{u))C[ip] dxdt. 


Since \\un\\L°°(Qj.) is uniformly bounded, ipn ^ T locally uniformly in K by assump¬ 
tion (ii), and |<p„(m„)| < |(/ 3 „(u„) — ip(un)\ + |<p(u„)|, we obtain for n sufficiently 
large 


(4.1) \\ipr,{un)\\L«>(QT) < sup{|(/3(r)| : |r| < C} -k 1 =: C^. 

Then, using assumption (i), we get 

^ Cip [ [ \C^"’[ip] —C[ip]\dxdt ^ d 

Jo Jr" 

as n —>■ oo. By the uniformly boundedness of and since ipn ^ T locally 

uniformly in M by assumption (ii), 

Ilv 5 n(u„) - (^(m„)||l~(q.j,) < sup{|(^„(r) - I^(r)| : |r| < C }0 as noo. 


(fniun) {C^^"^[ip] - /JW\) dxdt 
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Since we assume that £[ 1 /)] G L^[Qt), 


'0 


{^Pn{Un) - ip{Un))C[lp] dx At 


< \\Vn{Un) - (/j(Un)||L°°||£[l/']||Li 0 


as n ^ oo. By assumption (iii) and (A;^), \^p{un) — </5('w)| —>■ 0 a.e. in Qt as 


n —>■ (X), and ||‘/3('u„)||icxj(Q^) < C for some C independent of n. Hence, \^iun) — 
(/3(u)| is bounded by 2C. Moreover, since £[■)/)] G L^{Qt), Lebesgue’s dominated 
convergence theorem yields 


/ / {(fiiun) - (piu))C['tp]dxdt 

lo Jr" 

as n —>■ oo. The proof is complete. 


< / / |(^(u„) - (/3(it)||£['!/’]| dxdt0 

Jo Jr" 


□ 


Let us turn our attention to proving the other main results in this section. 


Theorem 4.1. Assume (A^), (f G ^3(0) = 0, and uo,uo G 


(a) There exists a unique entropy solution u G L°°{Qt) H (^([O, T]; of 

(O). 


(b) Ifu,u are entropy solutions of (1.1) with initial data uo,uq respectively, then 
for all t G [0, T] 

||M(-,t) - M(-,t)||Li(B«) < Iko - UoIIl1(R«)- 

(c) If u is a entropy solution of ([TtJ with initial data uq, then for all t G [0,T] 

lk(•)^)llz,l(R”) ^ I|wo||li(rn) and ||m(-, t)||ioo(R]v) < ||uo||loo(r]v). 

Entropy solutions are defined in Definition 2.1 in [23], and the result holds by 
Theorem 5.5 in [23] and Theorem 5.2 in [22] • 

In what follows, we let uq G n and define 


(4.2) ‘Prj{x) := (p * uJriix) — ip * oJrjiO) where w,, is given by (|1.7|) with TV = 1. 


Hence tpr, G W^ioc°°( 


C C'(M), it is nondecreasing by (A,^), (pniff) = 0, and 
iprj ^ T locally uniformly in M. Let be the entropy solution of (1.1) with 
replacing ip. Since entropy solutions are distributional solutions (cf. Theorem 2.5 
ii) and Section 5 in [22]), 

[ [ (urjdtip+(prj(ujj)C^[ip]) dxdt+ [ Uo-iplt^o dx = 0 Vi/'G C'“(M^x[0,T)). 

Jo Jr" ^ ^ Jr" 


Going to the limit as ly —> O’*" in (4.3), we will prove the existence and the a priori 


results given in Theorems 2.8 and 2.10 


Remark 4.2. We will prove that the L^-contraction holds for limits of the functions 
{ur;}r;>o- As a consequcncc of uniqueness (Corollary |2.4| ), this result then holds for 
all n L^-distributional solutions of 

Before these results can be proven, we need an auxiliary lemma. 

n £1 


Lemma 4.3. Assume (A,j), uq G L' 


, satisfy (A,^) for all 


r] > 0, and iprj ^ ip locally uniformly as ry —>■ 0+. If Urj solves (4.3) and satisfies 
Theorem 4-1 (b) and (c), then there exists a subsequence and a u G 

(^([O, T]; L^(IR^)) such that as rjn —> 0"^ 

Ur,„^u in C'([0,r];£j(,^(M^)). 
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Moreover, for all t S [0, T] 

lk(•)^)II li(r«) < I|uo||li(rn) and ||M(-,t)||2,cx=(R]v) ^ ir«ollL“(R»)- 

Proof. We will use Kolmogorov’s compactness theorem in the form of Theorem A .8 
in [5S]. Let K C be any compact set. 

(c). 


4.1 


Step 1: Urj is bounded independently of rj in Qt by Theorem 
Step 2: Since O) is translation invariant, v{x,t) = Uri{x + h,t) solves (|4.3[) with 


initial data Vq{x) = Uq{x + h) for every h S Let 7 S K". By Theorem 
and since translations are continuous in , 

sup / \ur^{x + h,t) — Urf{x,t)\dx < sup / \ur^{x + h,t) — Urf{x,f)\dx 
\h\<\'i\JK |/i|<|7|-/r~ 

< sup / \uo{x + h) - uo{x)\dx < max A„o(|/i|) =: 

|/t|<|7|jR» I'»l<l7l 

for some moduli of continuity A„„, A„p. 


(b) 


Step 3: Let ojs be defined by (1.7) and let 0 G C'“((0,T)). For any x G take 
ip{y,t) = <d{t)ujs{x — y) as a test function in (4.3) to find that 

rT 


0 = 


(4.4) 


iur,{y, t)ujs{x - y)Q'it) + (pr,{ur,{y, t))£^[u}s]ix - y) 0 (t)) dy dt 
({uni-, t) * ujs){x)Q'{t) + t)) * C^[ws])ix)e{t)j dt. 


lo 

pT 


For ps{t) defined by (|1.8[), we choose 


0(t) := <dg{t) = f (pg{T - ti) - pg{T - ta)) dr, 

J —00 ^ ' 

where 0 < ti < ^2 < T’- For 5 > 0 small enough, 0^(t) is supported in [0, T] and is a 
smooth approximation to a square pulse which is one in and zero otherwise. 

By ([ 43 , 

/ Ps{t-t2){urf{-,t)*0Js){x)dt= / Pg{t-ti){urf{-,t)*UJs){x)dt 

Jo Jo 


Let u°^{x,t) := Ujj{-,t) *ujs{x). By Theorem 4.1 (c) and the properties of mollifiers, 
we send 5 —>■ O’*" in the previous equality to obtain the following pointwise identity, 

(4.5) ufi{x,t 2 )-ufj(x,ti)= f {iprj{Ur,i-,t)) * Cf"[uJs]){x)dt. 

Jti 

Now, we need to estimate the integral involving the mollified version of Let 
t, s G [0, T] and take 6 < min{t, s}. Use ( |4.5| ) to find that 

]{x,t) - u^^{x,s)\dx < [ [ |((/?^('u^(-,t)) */:'^[a; 5 ])(a:)| drdx 
Jk Js 

= f f f IPrii^rii^ - y,T))\\C^[oJs]{y)\ dydxdr 
Js JkJw" 

< \\iPr,{Ur,)\\L<^(QT)\\^^[aJs]\\L^m)\K\\t - s|, 

where |iG| denotes the Lebesgue measure of the compact set K. As in the proof of 
Theorem |2.6| (see ( |4.1[ )), we obtain for ry sufficiently small 

\\Pviuv)\\L--iQT) < sup{|(^(r)| : |r| < ||uo||l~(r»)} + 1. 


> K 
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Moreover, Lemma 3.5 (b) and the properties of mollifiers yield 


ll>C'"[a; 5 ]||ii(Rjv) < (d + 1) / min{|zp, 1} d/r(z). 

J |2;|>0 

Hence, taking <5^ := |t — s|3 we see that 

(4.6) J lu^ix, t) - ufj{x, s)\dx <Ck, ip,no,p, {\t - s\^+ \t-s\^ , 

where 

CK,ip,UQ ,/j, 

= \v>{r)\ + l) [ min{|0p,l}dAi(z). 

^ ^ ^ IrKlluollr.oo ^ j\z \>0 


By the triangle inequality and Theorem 4.1 (b), 

/ |M,;(a:, t) — s)| da; 

Jk 

< / |u,,(a;, t) — u^(x, t)| da; + / |u^(a;, t) — M^(a;, s)| da; 

Jk Jk 

+ / |u^(a;, s) — u^(a;, s)| da; 

J K 

< sup t) - M^(-+ cr, t)llii(Rjv) + / |■u^(a;,^)-■u^(a;,s)|dx 

|cr|<5 Jk 

+ sup ||u^(-,s)-M^(-+(T,s)||il(RN) 

|cr|<(5 

< 2 sup ||uo - uo(-+ cr)||Li(R«) + / |u^(a;,t)-u^(a;,s)|da; 

|<t|<5 Jk 


< 2 max A„q (<5) + 


, ' / \u^ix,t)-u {x,s)\dx, 

\<^\<^ JK 


where is defined in Step 2 . Hence, by (4.6l 

J \urj{x,t) - u^(a;,s)|da; < A„o (|t- s|3^ +Ck, p,uo,p (l^- s|3 + |t- s|) 

— ■ J^K,p,uo,pi\t ~ s|) 

for some moduli of continuity A„(, and Ak,p,uq,p- 

Step 4: The assumptions of Theorem A.8 in [29] hold by Steps 1-3, so we conclude 
that there is a subsequence such that 

Ur,„^u in C'([0,T];L(„,(K'^)) 


as 77„ —>■ 0“*’. Finally, u inherits the properties of given in Theorem 4.1 (c) by 
Fatou’s lemma, and the fact that the limit of a uniformly bounded sequence which 
converges a.e. is also bounded. □ 

Remark 4.4. If was not fixed in the above result, but rather /a = (with /a„ 
satisfying (A^)), then the result still holds and the proof is the same provided we 
also assume that for some a > 0 there exists a function / S T“j,((0,a)) such that 

ll'C'"’'[w5]||ii(RN) </((5) for every de(0,a), 

where ojs is defined by O. Observe that the above ineq uality follows from the 
assumption sup„min{|zp, 1} d/a„(z) < oo in Theorem 


2.12 


Now, the proofs of the existence and the a priori results follow. 
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Proof of Theorem \2.^ Let be the solutions of (4.2) (cf. Theorem 4.1), u G 
L°°{Qt) n L^{Qt) n C([0, T]; the function provided by Lemma 4.3 and 

define := (that is, := = £), := and u„ := Ur,^. Then 

assumptions (i), (ii), and (hi) in Theorem 2.6 are satisfied by the n-independence 
of (|4.2|, and Lemma |4.3| Moreo ver, s up„ \\un\\L<^j^) < ||uo||l“(r«) < oo by 


Theorem 4.1 (c). Hence 


3 y Theorem 2.6 u satisfies (1.1) in the sense of distribu¬ 


tions: cf. Lemma 2.21 and Definition 2.2 Moreover, we have that u — uq G L^(Qt)- 


So, u is in fact a distributional solution of (1.1) in the sense of Definition 2.2 and 
it is unique by Corollary |2.4| 


Thus, any subsequence has the same limit, and hence, the whole sequence 
{uti}ti>o converges since it is bounded by Theorem 


4.1 


(c). 


□ 


Proof of Theorem 2.1(][ (a) Let u,, be the entropy solution of (4.2) (cf. Theorem 
4.1). Using the semi entropy-entropy flux pairs 


(ujj - ky 


and ± sign^(u^ — k){f{un) — f{k)) for all k G 


and the corresponding definitions for entropy solutions in stead of the Kruzkov 
entropy-entropy flux pairs in [22], we obtain 

/ {uTj{x,t) — dx < / (mo(x) — mo(3 ^))~''dx 

Jr" Jr" 

for Un,Urj G L°°{Qt) C L^{Qt) H C([0, Tj; with initial data uq^uq G 

n See [25] for the result and a proof. 

By Lemma [4.3[ we can take subsequences such that -G u,u a.e. in Qt 

as > O'*'. Thus, Fatou’s lemma yield the result. 

(b) By the contraction estimate obtained in part (a) and mq < Uq 8 ..e. in for 
all t G (0, r), Jgjv {u{x, t) — u{x, t))+ dx < 0. Hence, {u — m)+ = 0 and u < u a.e. 
in Qt- 


(c) Follows by Lemma 4.3 


(d) Follows by Lemma 4.3 


(e) Using the triangle inequality, and taking u,Uri„ as in Lemma 4.3 we obtain by 


Step 3 in the proof of that lemma that for all t, s G [0, T] and any compact set 
K 

||u(-,t) - u(-,s)||il(K) 

< \\u{-,t) - M^„(-,t)||il(^) -k - Ur,^{-,s)\\Li(K) + I|m,,„(-,s) - u{-,s)\\l^K) 

— 2||m(-, i) — t)||c([o,T];Lf^^(K«)) + J^K,cp,uoTi\t ~ sl) 

for the modulus of continuity (see the above mentioned proof). Since 


-G u in C^O, T]; Ll^^fR^)) by Lemma 4.3 the proof is complete. 


(f) Consider a standard cut-off function 0 < X G C^{R^) such that X{x) = 1 
for |a;| < 1 and X{x) = 0 for |x| > 2. We will write Xji(x) = T’(^) for R > 0. 


Following the proof of Lemma 3.10 (b), with 8a as defined there, we can take 


if{xR) = Xfi{x)9a{t) for any i? > 0 as a test function in Definition 2.2 (cf. Lemma 

Oait) 


2.21). Hence 

1 r 


/ u{x,t)Xii{x) dxdt = / 

a Js-a Jr" Jo 


-k 


IR" 


ip{u{x,t))C^[Xii](x) dx dt 
uq{x)Xr{x) dx. 
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Since Xu is compactly supported and u S C'([0, T]\ we can pass to the 
limit as a — >■ O’*' in the first integral to get /jjn u{x, s)Xii{x) dx. For the second in¬ 
tegral, we know that (p{u) G L°°{Qt), C^[Xii] G and 9a —>■ l[o,s) pointwise 

a.e. as a —>■ 0+, and thus, it converges as a —>■ 0+ to V’iuix, t))C^[Xp(\{x) dxdt 

by Lebesgue’s dominated convergence theorem. In this way, we get 


f u{x,s)Xji{x)dx = f f (p{u{x,t))C^[Xii]{x) dx dt + f uq{x)Xpi{x) dx. 

ilRW JO JR" JR" 

The function Xr converges pointwise as i? —)■ oo to 1, and it is also bounded 
by 1. Then, since u{-,s),uo G Lebesgue’s dominated convergence theorem 

allows us to pass to the limit as i? —>■ oo in the first and the last integrals to 
get u(x, s) dx and fgjv uo(x) dx, respectively, for all s G (0,T). Consider the 
nonsingular part of the Levy operator, i.e., (x + z) — X^ (x) dfi(z) which 

is bounded by 2/j,({z G : \z\ > 1}) for every x G Since Xfi{y) —>• 1 pointwise 
as i? —>■ oo for all y G Lebesgue’s dominated convergence theorem shows the 
pointwise convergence to 0 of the nonsingular part. For the singular part. Lemma 
3.5 (b) gives 


'0<|z|<l 


Xr {x + z)- Xr (x) dy(z) 






'bl<i 


■dy{z) 


which also goes to 0 as i? —)■ oo. Moreover, by the assumption |<p(r)| < Ls\r\ for 

kl < s, 


\\‘P{u{x,t))\\LHQT) < 


'|u|<5 


Ls\u{x,t)\dxdt+ ||(/j(m)||l=o(q,j 


dx dt. 


' >5 


Since u G L^{Qt), both terms on the right-hand side of the estimate above are 
finite. Then by Lebesgue’s dominated convergence theorem. 


[ f (p{u{x,t))C^[XR]{x)dxdt 
Jo Jr" 

The proof is complete. 


as 


R 


00 . 


□ 


5. Applications of stability 

This section focuses on proving the results stated in Sections [2.1| and [2?^ 

5.1. Compactness, local limits and continuous dependence. 

Proof of Theorem \2.1^ (a) Note that the sequence of solutions {ura}„GN satisfy the 
hypothesis of Theorem |2.10| By the assumptions. Remark |4.4[ and Lemma [T3l the 
result follows. 


(b) This is a consequence of the stability given in Theorem 2.6 For the initial 
condition, note that by the assumption sup„ ||wo,n||L<=°(R«) < oo and Fatou’s lemma, 
Uq G n and the convergence of Uo^n{x)'4’{x,d) dx follows by 

the L[Qj,-convergence of {wo.njnGN- D 


Lemma 5.1. Assume (), s G (0, 2), = —(—A) 2 , and if G C[ 


'OO rm,N\ 


Then 


lim \\ — {—A)2ip — Aif\ 


s —^2 


Li(E^) 


= 0 . 


Proof. The fractional Laplacian has a representation in the form (1.31 with measure 

-1 


dy = CN,s 


dz 


\N+s 


for 


cn,s = N 


1 — cos(zi) 


\N+s 


dz 
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where 
(5.1) 

see e.g. Proposition 4.1 in [33]. Hence 


-{-A)i^l){x) = CN, 

+ C_/V 


lim Cat s = 0, 

s —^2 ~ 


ip{x + z) — ipi^) — z ■ Diplx) 


|N+S 


dx 


ld<i 1^1 

f tj;{x + z) — ip{x) 

\z\>l 


|N+s 


dz, 


where the last term goes to zero in (M.^) as s —>■ 2 since it is bounded in 
by cjv,s2||V'||li(r«) dz for s > 1. 

The explicit form of cn^s given in (5.1) yields 

Atp{x) = A'iI){x)cn,sN ( I dz + AV’(x)CAr,siV ^ ^ cos(zi) 


'|z|<l 


ld>i 


|Af+s 


d0. 


Again, the last term goes to zero in as s —>■ 2 since |1 — cos( 2 :i)| < 2 and 

then it is bounded in by cn,s‘ 2N dz for s > 1. Using Taylor’s 

theorem, we see that 1 — cos(zi) = ^zf — ^zf cos(^) for some ^ € [0, zij. Hence, 


, , f 1 — cos(zi) , 1 

/ |r|Af+s ~ 2 / 121-^+^ 

and the following estimate holds: 


“ oT cos(C) / 

j\z\<l 


|N+s 


dz, 


N f 

— AiI}{x)cn,s / 

J\z\<l 


cos{^)i 


I A+s 


^ dz 


Li(R«) 


N f 

< ^Cjv,s||AV’||li(m«) / 

J h 


|z|<l \Z 


N-2 


dz 


which goes to zero since dz < oo and (5.1) hold. 

To estimate the remaining term in (5.2), note that for all r > 0, 


and then 


NAtp{x) [ z^ dz = A^(a;) f |zpdz= f D‘^'tp{x)z ■ zdz, 
J\z\<r J\z\<r J\z\<r 

^D‘^jp{x)z ■ z 


^CN,sNAtp{x) 


^ dz = CN,s 


i\z\<i i^r+ 

We combine the all above estimates to get 
lim 11 - (-A)5'!/'(a:) - AV'(a:)||ii(R«) 

S-S.2- ' ' 


'ld<i 


I A+s 


dz. 


= lim Cat , 

s — h2 ~ 


iIj{x + z) — ^l;{x) — z ■ D'iIj{x) — ^D‘^iIj{x)z • z 

^-dz 


'ld<i 


< lim_ CAr,s-||T>^t/’llLi(»") 




lA+s 


dz, 


Li(R») 


where the last inequality follows from Taylor’s and Fubini’s theorems. Since the 
z-integral is bounded by dz < oo and (5.1) hold, the limit is zero and 

the proof is complete. □ 

Proof of Corollary \2.13[ (a) We will use Theorem |2.12 and R ema rk 4.4 to prove 
the result, and now we verify the assumptions. By Lemma 5.1 —(— A)i-0 —>■ 
Alp in as s ^ 2~ for all ip G Moreover, by Lemma 3.5 (b). 
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properties of mollifiers, lims_^ 2 - /| 2 |<i 
(see previous proof), 


|2 Civ.s dz 


= 1, and lim,^ 2 - 4|>i fjyw'' 


{-A)iu:s\\L^ = hlD^LOsh^ f 

J h 




<C[1 + 


1 


ld<i 


p|«+ 

I ^|2 ^A,s d2: 
A+s 


= 0 


■ 21^511 


Li 


'|z|>l 


cn,b dz 

\N+s 


for s clos e to 2. Henc e, si nce also ip is fixed (independently of s), we may use Theo- 
1 and Remark|4.4|to get a subsequence and au G (^^([0, T]; Ljq^(K^)) 


rem 


2.12 


such that Usj —?► rt in (^([O, T]; Ljqj,(M^)) as j —?> oo. Finally, the uniqueness for the 
limit (equation) [18], and the boundedness of the sequence {us}sg(o, 2 ) (Theorem 


2.10 (d)), ensures that the whole sequence converges. 

(b) Since (—A)t'0 — {—A)iil; in L^(K^) as n —>■ oo (a similar argument as in 
Lemma 5.1), = r"*" —>■ Pfn{f) = u™ locally uniformly as n —>■ oo, and 

|j(—A)T^a; 5 ||ii < (7(1 + 5~'^) by the proof of part (a), convergence for a subse¬ 
quence follows by Theorem |2.12[ Moreover, the convergence of the whole sequence 
follows from uniquenes of the limit (Corollary |2.4| ) and boundedness of the sequence 
(Theorem 2.10| (d)). □ 

5.2. Numerical approximation, convergence and existence. We start by 
showing that a standard finite difference approximations of the Laplacian can be 
written in the from (1.31 and that convergence of the resulting scheme then follows 
from our theory. 

Example 5.1. Let G K" for i = 1,...,A be points with i-th component 1 and 
the other components 0. Using J-measures and > 0, we define 

N 




2=1 


^hei + 5- 


hci 


It is clear that pt is a measure satisfying ( |A^| for every h > 0. Moreover, 

N 


C^''[v]{x) := f v{x + z) - v{x) dfj,h{z) ='^ 


v(x + hci) -|- v{x — hci) — 2v{x) 

W 


can be reformulated as 

,,{x ■ 


With p = ph, problem 

/.ON n / ./N ^piuh{x + hei,t)) + p{uhix-he^,t))-2p{uhix,t)) 

(5.3) OtUh{x, t)- 2_^ -- = 0 

in V'{Qt). 

For Ip G an application of Taylor’s theorem reveals that there is a 

C > 0 such that 

f \C^''[ip]{x) — A'tlj{x)\ dx < h'^C\\D'^'4>\\i^i(^s^N-^ ^ 0 as /i —>■ O'*". 
Moreover, for h small enough, 

sup / mm{|z| , 1} dph[z) = sup > - ^ - = 2N 

h J\z\>o h j-pp h 

Then by Theorem 


i=l 


2.12 


there existis a subsequence {uhj}jen of solutions of (5.3), 
and a u G (7([0, T]-T^(IR^)) such that tt/i, —» u in (7([0, T]; L[qj.(K^)) as j —>■ oo. 
Moreover, the limit u satisfies equation (1.5): 

dtu — Ap{u) =0 in V'^Qt)- 

In fact, as in the proof of Corollary 2.13 the whole sequence {uh}h>o converges. 
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We can proceed as in this example to get convergence for a more general class 
of second order local operators. 

Lemma 5.2. Assume /i > 0, P G N, u = (cti, up), Ui G for i = 1,P, Lf, 
is defined by (2.9), and if G C“(]R^). Then 

LlW{x)= [ U{x + z)-'if{x))d^j,h{z)=:C^^''’‘’[if]{x), 

J\z\>0 ^ ' 


where the measure ^ P^-ha-i)- Moreover, satisfies (A^ 


sup/ min{|zp, 1} d^h.cr(z) < oo, 
h J\z\>0 


and 




as 


] - tr[crcr D ifWlpi 0 
Proof. By an elementary identity and Talyor’s theorem, 

p 

D^ifi^x)] = '^{ui ■ D)^'fi{x) 


h^0+. 


= E 


2=1 


v{x + hai) + v{x — hai) — 2v{x) 

■ 


N 




/3! 


i=l |/3|=4 ' 

Here we use standard multiindex notation, with multiindex /? = (/3i,..., fin) G N'^, 
to account for the 4-th order derivatives. Since the first term of the last line is 
Lf[ijj]{x), the rest of the proof follows along the arguments of Example 5.1 □ 


We aim to consider the general operator defined in (1.3). In order to use our 
stability result, we would like to prove that the operator Lf, defined in (2.10) is a 
particular case of the operators studied in this paper. The following result ensures 
this fact. 


Lemma 5.3. Assume (A^), h>d, Cl is defined in ( |2.10[ ), and if G C° 
Then 

cimx) = 


|z|>0 


x + z) - Ip{x)'^ dvh{z) =: ['fi]{x) 


where the measure Vh = Sa/o T (zid^ + Rh) ■ Moreover, Vh satisfies (A^ ) and 
sup / min{|zp, 1}dz//j(z) < oo. 

h J\z\ >0 

Proof. By the definition of (5^^, it immetiatly follows that C^ = C"'^. It remains to 


show that Of, satisfies ( |A^p . For h < l/\/N, 

J ^dPh{z)= t^{Za+Rh)<tJ-(J^\z\>l-VN^yj <fJ.(J^\z\>^'^ 


b„|>l 


which is finite since fi satisfies (A^). Moreover, for h > 0 small enough, 
|zp dvhiz) 


J\z\<l 

^ E 


0<|zq|<1 


dt^iz) < E / 

0<|z„|<l''^“+''^'* 


< 


'?i/2<|z|<l + vW| 


z\+VN^'^\yciz) 

zI + Vn'^'] d/i(z) < (l + Viv)Y \z\^dyiiz), 

/ »/ 2 2 


which is also hnite since y. satishes (A^). The proof is complete. 


□ 
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Lemma 5.4. Assume (A^), and are defined in (1.3) and (2.10) respectively, 
and ijj G Then 


^ 0 as /i^0+. 

Proof. The following inequality is just a use of the definitions, 


/b» 


l^hbPKx) - ^^bP]{x)\ dx 


dx 


< 


[ '^{ij{x + z^)-fi{x))[ dp,{z) 

Za-\-Rh 

(^{x + z) - fi{x)^ dpL{z) 

X + z) — dp{z) 

X + Za) — ifix + z)^ dfj.{z) 


a/0 

aGZ” 


E 

ct#0' 


' Za+Rh 


dx. 


We will show that both terms go to zero with h. Indeed, for \z\ < 1 we have that 
\z\^\ii,^{z) -G 0 pointwise as h — » 0+. Then, by Lebesgue’s dominated convergence 
theorem, (), and Lemma 3.5 (b), we have as h —>■ O'*' 


[ [ ('^{x + z)-'ilj{x))dfi{z) dx<h\D^ filial [ |zplj^^(z)d/i(z) ^ 0. 

Jr^ jRh ^ J I 9bl<i 

For the second term, we need to consider separately the cases when when we 
are close or far from the origin. First note that for any z G Zad- Rh we have that 
\za — z\ < VN^. Since p, satisfies (A^) and if G 


Avt *— 


/R« 


/ (ljj{x + Za) - ifix + Z)] dfi{z) 

\a\h>l J^c+Rh 


dx 


< I1-Di/’||li(h«) / \za-z\dp.{z) 


< ^^ll-OV'IUbR") / d/r(z)0 as /n> 0+. 

On the other hand, by the symmetry of /r and also of the term in the sum, we have 
that 



dx 

z) ■ Dip{x)'^ 


dfi{z) 


dx. 


We make use of the Taylor expansions 

lf{x + Za) = lf{x + Z) + Dlp^X + z) ■ (Za - z) + Gi(x,Z, Za){Za “ z) ■ {Za - z) 
DiIj{x + z) = Dif{x) + G2(x, z) ■ z 
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where ||G'i||ii(RN^da:) + l|G 2 ||Li(R«_d:c) < C\\D'^tf}\\Li(^N) for some constant C > 0. 
In this way, 


Iint< / X! / (G2Z-{Za-z)+Gi{Za-z)-{Za-z))dfi{z) 


da; 


f (I 

zjjZa - -Z| + jZa “ 

Za~\~Rh 

/ 


d^i{z) 


< G\\D'^'4’\\l^m") [ + d^(z) —> 0 as h—>-0^. 

a |<|z|<i+%/]v| 2 \ 2y 

Since the integrand is dominated by 2|zp which is an integrable function with 
respect to the measure ^ on the set {z £ : \z\ < 1} by (), the last term goes 

to zero by Lebesgue’s dominated convergence theorem. □ 


Proof of Proposition \2.16\ Note that by Lemmas 5.2 and 5.3 and are in the 


class of operators defined by (1.3) and (A 


(a) Existence, uniqueness and regularity follow from Theorem 2.8 


(b) Follows from Theorem 2.10 (c) and (d) and interpolation. 


(c) Lemmas 5.2 and 5.4 ensure the -consistency. 


(d) Follows from Theorem 2.10 (b). 


(e) Follows from Theorem 2.10 (f). 


□ 


Proof of Proposition By Lemmas 5.2 and 5.3 and Proposition 2.16 

sup/ miTL{\z\'^,l}d{pLh,a + Vh){z) <oo, 
h J|2|>0 


.cj, 


and 


(A 




+ [■!/;] II —>■ 0 as 


O’* 


Since also and ug are fixed (that is, independent of h), by Theorem 2.12| there is 
a subsequence {M?i„}nGN of solutions of (2.5), that converge in C([0, T]; 
to a function u. Moreover, this function u is a distributional solution of ( |2.6[ ). 
Finally, u also belongs to L°°{Qt) H L^{Qt) by Proposition |2.16| (b) and Fatou’s 
lemma. □ 


Proof of Corollary \2.18{ A ny limit point u from Proposition |2.17] is a distributional 
solution of (2.7) and (1.2). □ 


Proof of Theorem By Proposition |2.17 there is a converging subsequence with 
a limit u which has the right regularity and is a distributional solution of (2.6). As¬ 


sume there is a subsequence that converge to another limit v. Then by Proposition 
|2.17| again, there is a subsubsequence that converge to a limit which is a distribu¬ 
tional solution. By uniqueness of the limit, z; is a distributional solution. But then 
V = uhy the uniqueness given in Corollary |2.4| for the case ct = 0 or the local result 
in [18] . Hence all subsequence limits are equal to u and since the sequence itself is 
bounded (Proposition 2.16 (b)), the whole sequence converges to it. □ 
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6. Auxiliary elliptic equation 


In this section we study the elliptic equation (3.1) introduced in Section with 
the ultimate goal to prove Theorem |3.1| and Lemma |3.4| We will also need the 


following approximation of (3.1) where the measure /i is replaced by := 


( 6 . 1 ) 

with e > 0, 


evs^rix) - C^^[vs^rKx) = g{x) in 


cN 


[iIj]{x) = [ (il}{x + z) -'tp{x)) dgr{z). 
J\z\>0 ^ ' 


Note that for any r > 0, the operator £^'■[-0] is well-defined for merely bounded 
Ip, and that Lemma [3.5| also holds for see Remark 3.6 (b). Also recall the 
notation = {el — C^)~^ and define BP’- := {el — £''’’) 


Remark 6.1. (a) Since (3.1) and (6.1) are linear equations, we have formally, for 
any multiindex a € NW that is a solution of (3.1) or (6.11 with right hand 


side D°‘g if u is a solution of the same equation with right hand side g. 
(b) Let Ip S C'^(K^), and let p G {1, oo}. Since 

{C^ - C^'^)['ip]{x) = [ (^^{x + z) -'4){x) - z ■ Dpj{x)] dp,{z), 

J\z\<r ^ ' 


we have that £'^'■['0] —)■ C^[ip] in as r —>■ O’t by Lemma 3.5 (b) and 

Lebesgue’s dominated convergence theorem. 


6.1. Preliminary results. We will state and prove a very general Stroock-Varopoulos 
type of inequality which is of independent interest. First we consider the bounded 
operators £^''. 


Lemma 6.2. Assume ( A^ ), tp G L°°{M.^)nL^{M.^) and C G (^(K) is nondecreasing. 


Then for any r > 0 we have, 


Ir ■■= [ C{'tp{x))Cf^’-[lp]{x)dx 

= -\[ [ {C{tp{x + z))-C{'^p{x))){^p{x + z)-'^p{x))d|J.r{z)dx, 

^ Js." J\z\>0 


and in particular, A < 0. 


Remark 6.3. More generally, the above lemma holds as long as the integral A is 
well-defined for ip and C(V')- 

In the proof we need a technical lemma which wil be proven in Appendix fA} 

Lemma 6.4. Assume n is a nonnegative, symmetric and locally finite Borel mea¬ 
sure on . Let A,B be Borel sets on R^, and let 


Mi{A,B) 


[ ( [ du(a:)] dz = f v{B — z)dz. 

JA \Jb-z J Ja 


M 2 {A,B)= f ( [ du(x)) dz = f iy{A — z)dz. 

JB \Ja-z j Jb 

Then Mi{A,B) = M 2 {A,B). 

Proof of Lemma [£l| Observe that (p{ip) G L°°{R^), and since |£^''['!/^]| dx < 
2 ||'*/'IIli /|z|>r. d/r( 2 ), £^’'[ip] G L^{R^). Hence A is well-defined. 
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By the symmetry of /i, the gradient term in the nonlocal operator vanishes. 
Fubini’s theorem and a relabelling of the variables gives 


Ir = 


/ C(V'(a;)) / (ipix + z) --ipix)) dnriz)dx 

/b« J\z\>0 ^ ' 

= / -'llj{x))ll^_^lyrdfj,-xiz)dx 

JR" i|z-x|>0 

= / / C{'>Piz)){tjj{x)-^p{z))llx-z\>rd^J,-z{x)dz. 

Jr" J\x-z\>o 

Since /i is a nonnegative, symmetric and finite Radon measure on (and hence 
a Borel measures), we can use Lemmato see that 


’ \x—z\>0 


C{tfj{z)){ilj{x) - i;{z))lix-z\>r dn-zix) dz 


Ir" -'|a;-z|>0 
It then follows that 


2Ir = - 


Ir" J |z-x|>o 


C{i^{z)){ll^ix) - 'tfj{z))llx-z\>r dfJ.-x{z) dx. 
C{ilj{x)){tfj{x)) - ilj{z))li^_x\>r dfi-xiz) dx 


’ lx —2| >0 


' I 2 :—x|>0 


C{iIj{z)){iIj{x) - 'ip{z))lix-z\>r dn-xiz) da; 

I 

(C(^/>(a;)) - Ci-fPiz))) (i’ix) - '<Piz))i\z-x\>r dfi-xiz) da;. 


Since (C('0(a;)) — ({tjj{z))) — V’(-z)) > 0 for all x,z G , Ir < 0. 

Now we give the general result, considering the general nonlocal operator 


□ 


Corollary 6.5 (General Stroock-Varopoulos). Assume (A^), and ( G C'^(K) such 
that > 0. 

(a) Letip Then 


I := 


(p{ip{x))C^[ip]{x) da; 


2 Jr" J\z-x\>o 


(Ci.i’iz)) - Cii’ix))) (V’W - ip{x)) dfj,{z) da; 


< 0 . 


(b) Let Ip G n If Z G C'2(M) is such that Z{0) = 0 and 

[Z'Y = C,', then 

[ C{'>Pix))^^[tp]{x)dx ^ [ [ {z{'ip{z)) - Z{'ip{x))Ydn-x{z)dx. 

Jr" ^ aR« a|2:-a;|>0 

Moreover, 

(z{iP),C^^[Zp<P)]\ =-]- [ f {Z{Pj{z)) - Z{'tP{x))ydfi-x{z)dx 

\ / L^{R") ^ J|z-£c|>0 

= - {cnHz {^)]' . 


Remark 6.6. The (energy) norm in part (b) is much studied when CA" = —(—A) 2 , 
s G (0,2), and Z = I (see [7l[33]). In this case 


(fP 


(-A) 



L2(B«) 


1 

2 


(V’(z) - tpix)y 


' |2 —x|>0 


— -7-1 AT+s 


dz da; = 


(—A) i Ip 


2 

L 2 (RN)' 
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This is called the Gagliardo (semi)norm of ijj and is denoted by ,2 




Proof, (a) By Remark 6.1 (b) 


da; 


/r« 


< 




0+, 


and we may send r —>■ 0“*" in Lemma |6.2| to get 
['tf]{x) da; 


Jr" 

= — - lim 

2 r—>-0+ . 


[ (C(V'(2)) - C(V'(a;))) {ip{z) - tp{x)) d/i_^(z) dx. 

«/ I 2 —x|>r 


By the assumptions on ^, 1 /; and ( |A^[ ), (^(^/^(z)) — ^(i/;(x))) ('!/;(z) — > 0 is 

integrable with respect to d^-x{z) dx on x IR^\{0} since 


/r« 


' \z—x\>{) 


ici'fpiz)) - Citfix))) {tf{z) - ip{x)) d^i-x{z) 


dx 


( 6 . 2 ) 


J\z\<l 

+ 4||C(i/’)||l~||iAIUi [ 

J\z\>l 

Thus, Lebesgue’s dominated convergence theorem gives the desired result. 

(b) For a,b G R, the Fundamental Theorem of Calculus and Jensen’s inequality 
gives the following pointwise inequality: 

2 


(6.3) 


{Z{b)-Z{a)f = (^J^ Z'{t)d?j <{b-a)j^ {Z'{t)f dt 
= {b-a) f C{t) dt = {b- o)(C(6) - C(a)). 

J a 


By the assumptions, we can easily check that Z{'if) G 1F^’'^(M^) n 
So the integral 

If r 

2 Jr" J\z-x\>o 


{Z{tp{z)) - Z{tp{x))ydn-x{z) dx 


is well-defined using a similar argument as in (6.2). Then, part (a) and (6.3) gives 
the first result of part (b). 

Next, part (a) yields 

^ ^ ^ (Zitpiz)) - Z{ij;{x))ydfi_x{z)dx. 




L2(R«) 


IR" a|z-a:|>0 


Moreover, since Zlip) G kF^’°°(]R^) n then by Lemma 3.5 (b) and inter¬ 

polation, both Z{tp) and C^[Z{tjj)] are in L^(M^). We then conclude the proof by 
application of Lemma 3.7 and Remark |3.8| (b). □ 


6.2. Results for the approximate elliptic equation (6.1). We will now focus 
on proving some a priori, uniqueness, existence, and stability results for (6.1). 


Proposition 6.7. Assume (A^ 


(a) If g G L°°{R^) and Vg^r S L°°(R^) solves evg^r — Id^’'[ve,r] < 9 then 

^ll(■^e,r■)’'’||L“(R«) < II (fi')"''II l“ (R«) ■ 
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is an a.e. solution of (6.11, then 
') < ||5lU“(R«)- 


(b) If g € and Ve,r € L' 

(c) Let g,g,Vs^r,Vs,r G ev^^r - I^^"[ve,r] < g a.e. and eve^r-> g 

a.e. If g < g a.e., then Vg^r < Ve,r a.e. 

Proof, (a) Assume first that g,Ve,r G C'b(R^)- Then for all <5 > 0 there exists a 
xs G such that 

Ve,rAs) + > SUp{z;e,r}- 

Then, since is an a.e. solution, 

eVe^rixs) < gixs) + / (ve^rixs + z) - Ve,rixs)) dg.r{z) 
i|z|>0 ^ ^ 

< II(5)'^IU“(r«) + / ( sup{ne,r} - Ve,rixs)) dg{z) 

J\z\>r ^ ' 


"L ||(5)'''IIl“(r«) + g 




> ’’})• 


> ’’})), 


Hence, 

£:sup{w£,r} < £Ve,r{x)+£^ < II (ff)^ II L“(R«) + 5[£ + g{{z G : |Z| 
and we pass to the limit as d —)■ 0+ to get 

£sup{z;e,^} < ||(ff)+||L=o. 

In the general case, when g,Vg^r G we need a regularization argument. 

Let uf,. := ujs * Ve,r and mollify the inequality to see that 

£v^.,. — C^''[v^A^9s ia 

By the first part of the proof and the properties of mollifiers, 

'ye.r(a:^) < k£.r(a;)-?^£.r(a;)|+'yf,r(a;) < o(l) + i||( 5 )+||i=o(Riv) as d 0+ 
for a.e. x G M.^. Part (a) follows. 

(b) In a similar way as in (a), we find that esup{— < II (s)” IIl°°(r«) and 
combine with (a) to conclude that e||w£^r-||L°°(R™) < II<?IIl“(r«)- 

(c) Since w = — ^E.r solves ew — [w] < g — g, by (a) and the assumptions, 

it follows that £sup{w} < ||(g — g)''’||L°o(RN) =0. □ 


Proposition 6.8 (Existence and uniqueness). Assume (A^ 


(a) If g G C'b(R^), then there exists a unique classical solution r G Cb(R'^) of 

(O. 

(b) If g G then there exists a unique a.e. solution Vg r G of 

(1^. 

(c) If g G then there exists a unique a.e. solution G L^(R-^) of (6.11. 

(d) If g G then there exists a unique classical solution Ve,r G C^{RP)of 

(6.1). Moreover, 

sWD^^v.AIl^ < ll^“5lk- 

for each multiindex a G N'^. 

Proof. The proofs of (a), (b), and (c) follow from standard arguments using Ba¬ 
nach’s fixed point theorem. Let X denote any of one of the spaces C'b(K^), L°° (K^), 
and and note that A is a Banach space. Let the operator T be such that 

(6.1) is equivalent to the fixed point equation T[m] = u: 


T[v,]ix) := 


^ + I\z\>r 


/ Ve{x + z)dg{z) + g{x) ] . 

yJ|z|>r J 
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It is easy to check that T is a bounded linear operator on X, and straightforward 
computations also shows it is a contraction: 

\\T[ve]-T[v^]\\x <a\\ve-Vs\\x for a = ^ < 1. 

£ +jp|>r 

Hence by Banach’s fixed point theorem there exists a unique G X such that 
Ug = T[ug] in X and then also a.e. (everywhere ii X = Ci,)- 

(d) Let [ 5 ] and define Si^h'4’ by 

ijjix + hei) - iIj{x) 

Oi,hipw = --• 

By part (a), we have uniqueness for C'b(K'^) solutions of 
hci) = + hei)]{x), and then by uniqueness and linearity Si^hV = B^^[di^hg]- 

In addition, there exists a unique Wi^e,r G C'b(M^) such that Wi^e^r = B^''[dxig]- 
Using linearity and Proposition |6.7| (b), we get 

^i,h'^€,r\\L^ ^ ^i,hg\\L°^ - 

When /i —>■ O"*", Si^ng —>■ dxtg uniformly on and hence 6i^hVe,r —> Wi^e,r in L°°. 
This implies that dx^Ve^r = Wi^e,r- Moreover, by Proposition |6 . 7| (b), 

II II 7^00 — II IULeT II ^ II II ■ 

A similar argument shows that for each multiindex a G D°‘Ve^r = B^^[D°'g], 
and hence belongs to C'b(K^). □ 


(6.1). Hence, + 


Corollary 6.9. Assume 

LUK^). 


andgG Cb(K^). If{g)+ G L^iR^), then {B^-[g\)+ G 


Proof. Note that g G Cb(ffi^) implies that (g)+ e C'b(ffi'^). By Proposition 6.8 (a) 
and (c), and the assumption on {g)~^, we have that have that B^^[g] G C'b(K^'') and 
B^’'[{g)'^] G L^(M'^)n(7b(K^) are the uniq ue cl assical solutions of (6.1) with right- 
hand sides g, {g)~^, respectively. Proposition 6.7 (c) ensures that B^^ [(ff)^] > 0 since 


{g)'^ > 0. In the same way, we get B^’'[—{g) ] G C'b(M^) and B^^[—{g) ] < 0. 

Adding the equations for [((;)+] and B^'^[—{g)~], and noting that ( 5 )+ — 
{g)~ = 5 G C'b(R^), we get 

e(Bn(5)+] - BM( 5 )-]) - [bM( 5 )+] - B^[(g)-]] = g. 

It follows that B^^[g] = [((;)+] — B^’'[{g)~] by uniqueness. We conclude that 

0 < (B^[g])+ < BM(g)+], and thus, (BMg])+ G L^(R^). □ 


6.3. Results for the elliptic equation (3.1). Now, we state and prove compar¬ 
ison, uniqueness and existence results for classical solutions of ([ 33 . These results 
will be obtained from the corresponding results for (16.11) and limit procedures. 


Lemma 6.10 (Comparison). Assume ( |A^[ ), < 7,5 £ and v^, Ve G C'^(M^)n 

B°°(K^) are solutions of (3.1) with right-hand sides g,g respectively. If g < g a.e., 
then Ug < Ve in 


Proof. Note that w = — Ve solves ew — C^[w] < 0, and hence, also 

ew[w] < IKB'" -/:'"'-)MIIl“(r«)- 

By Proposition 6.7 (a), it then follows that 


£||(w)+||ioo(RiV) < IKB'")H II io 


')■ 


Assume for moment that w G Cb(ffi'^). Then by Remark 

IKB'" - 0 as : 


6.1 


(b), 

0+, 
















36 


F. DEL TESO, J. ENDAL, AND E. R. JAKOBSEN 


and we conclude that w < 0. 

The general case follows by mollification: ws = w * ws (cf. (HH)) satisfies ews — 
^^[ws] < 0 and hence by the first part of the proof and properties of mollifiers, 


w{x) < ws{x) + |r(;(a:) — ■u: 5 (a:)| < 0 + o(l) as <5 —> 0“'' 
for every x G . The proof is complete. 


□ 


Corollary 6.11 (Uniqueness). Assume (A^), and g G Then there is at 

most one classical solution Vr G C^{ 


nL° 


of ([ 33 . 


Proof. If g = g a.e., then Lemma |6. 10 


gives Ve = Ug in . 


Proposition 6.12 (Existence and Stability). Assume (A^l, g G C{ 


□ 


, e > 0. 


(a) There exists a unique classical solution B^[g] = Ve G of (|3.1[). 


(b) Any sequence {ue,r„}neN of solutions of (6.1) converges locally uniformly to 
Ve = B^[g] of part (a) as r„ 0+. 

Proof, (a) Let 0 < 0+ as n —>■ 00 , and let v„ := Vs^r„ G be the 

unique solution of (6.1) given by Proposition 6.8 (d). Moreover, for all n > 0, 


e|kn||L“ < IIsIIl”, e\\DVr,\\L^ < ||T>5|U-, 
epVlU- < WD^gh^, eWD^nh^ < WD^gh^- 

The sequences {wn}n> 0 ) {Dvn}n>o and {D^v„}n>o are thus equibounded and equi- 
lipschitz. By Arzela-Ascoli’s theorem there exists a subsequence (still denoted by 
Vm Dvn and D^Vn) such that {vn, Dvn, D^Vn) converges locally uniformly (and 
hence a.e.) as n —>■ 00 to a limit (v, Dv, D'^v) which is bounded and continuous. 

We check that Dv = Dv and D^v = D'^v. Let a G denote a multiindex. By 
Taylor’s theorem 


(6.4) 


yn(g) =v„(x) + Dvn(x) ■ (y - x) + -T>^v„(x)(y - x) ■ (y - x) 
+ - 7 ( 1 /- 2 ;)“ [ {f-t)‘^D°‘vnix + t{y-x))dt. 

Jo 


|a|=3 


Since 


^ —,iy - j (1 - t)‘^D°‘vn{x + t{y - x)) dt 

cd ./n 


|q|=3 


< -JD^gh^ E 


\y-x\ 


|a|=3 


al 


we can take the locally uniform limit in ( |6.4[ ) as n 00 to obtain that 

1 - 


viy) = v{x) + Dv{x) ■ {y - x) + ^D'^v{x){y - x) ■ {y - x) + o{\y - x\‘^) as y 
By definition, it then follows that Dv = Dv and D^v = D'^v. 


We now go to the limit in (6.1) as r„ —0+, and we may assume that r„ < 1. 


In order to show the convergence, the nonlocal operator in (6.1) will be written as 


[Vn]{x) = £ 1 ’'" [Vn]{x) + [ (vn(x A z) - V„(x)) dg(z), 
J\z\>l ^ ' 


Bi"" [Vn]{x) := / (vn{x + z) - v„{x) - z ■ Dv„{x)) dgr„iz). 
J\z\<l ^ ' 


with 
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By the triangle inequality and Lemma 3.5 (a), 

< - max + z) — + z)| f Izl^dnlz) 

2|ri<i' J\z\<i 

+ - max |Z?^i;(x + z)| f |zplu|<r d^{z). 

2 ril<i J\z\<i 


So, the local uniform convergence and Lebesgue’s dominated convergence theorem 
ensures that | [u„] (x) — £ 2 [^]( 2 ^)| —0 as r„ —0+ for all x G . The re¬ 

maining term in the nonlocal operator also converges by Lebesgue’s dominated 
convergence theorem: 


[ (vUx+z)-Vn{x)) dfi{z)f (vix+z)-v{x)) d^i{z) as 
a|z|>i ^ ^ J\z\>i ^ ' 


0 +. 
'Id; 

Sending r„ —>■ O’*" in ( |6.1[ ) then shows that v solves (3.1). Moreover, the limit is 
unique by Corollary |6.11[ 

(b) In fact, part (a) shows that all limit points of the sequences coin¬ 

cide by uniqueness (see Corollary 6.11). By Proposition | 6 . 8 | (d), every sequence is 


bounded, and hence, the whole sequence converge locally uniformly to the solution 
0 +. 


of (3.1) as T 


□ 


Proposition 6.13. Assume ( |A^[ ), g G c > 0, and = B'^lg], If 

{g)'^ G then 


: [ {vs)^ dx< f (g)+ dx. 


C^{R^) such that 


Proof. By Proposition 6.8 (d), for any r > 0, there exists a unique function G 


eVs^rix) - C^^[vs^r\{x) = g{x) in 
Consider fb G such that 0 < A, 

'l |x| < 1 


j,N 


Xix) = 


0 |x |>2 


and define Xr{x) = T’(|r) for i? > 0. Then for every r > 0, by Proposition 6.8 (d), 
there exists a function ur G (^“(IR^) such that 

(6.5) euR{x) — C^^[ur\{x) = g{x)XR{x) for all x G 

Let : K —>■ IR+ be a smooth approximation of the sign"*" function. More 
precisely, Cs{x) = 0 for x < 0, C,'g{x) > 0 and 0 < fs{x) < 1 f or x > 0. Since 
0 < (g<Ti^)+ < (g)+ G L^(IR''^), {ur)^ G by Corollary 6.9 and 


urCs{ur) dx 


gXRfsiuR) dx 


<ll(«R) + llLi||a(«R)l|L= 


< bill 


dx. 


|a;|<2R 


Then by ( |6.5[ ), [/^jv £'''' ['Rr]Ci 5 ('Wr) dx| < oo, and we may multiply ( |6.5[ ) by Q and 
integrate over to find that 


/r« 


URCs{uR)dx= / C^’'[uR]Cs{uR)dx+ gXRfs{uR)dx. 

Jr" Jr" 
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So, Lemma 6.2 and Remark |6.3| gives that < 0 and hence 

£ UR(s{uR)dx< / gXRCs{uR)dx< / dx. 

Jr"' Jr" 

Letting Csi^R) sign+(ufl;) as (5 —>■ 0+ in the above inequality (using Patou’s 
lemma on the left-hand side since urCs{ur) > 0 ) yields 

( 6 . 6 ) e f (uR)+dx< [ ( 5 )+dec. 

J'R" Jr" 

We note that the sequence {mr}r>o is equibounded and equilipschitz since 
Proposition 6.8 (d) gives 


1 , 


\\uRh^<-\\g\\L^ 
e 

WDurU^ < kn igXR) iu== < hogh^ + 

£ £ £\^i\ 

Hence, by Arzela-Ascoli, uh —>■ u as i? —>■ 00 locally uniformly in (and thus a.e. 
in K^). Sending i? —>■ 00 in ( |6.5[) shows that u = v^^r, that is, the unique solution 
of (3.11 given by Proposition 6.8 (d). Furthermore, we can send i? —>■ 00 in ( 6 . 6 ) 
(again using Patou’s lemma) to obtain 

e [ (^e.r)^ dx< [ ( 5 )+ dx. 

Jr" Jr" 

(b), we can let r„ — 0 + in the above 


6.12 


By Patou’s lemma and Proposition 
estimate to get 

e [ (u£)+ dx< [ (g)+ dx, 

Jr" Jr" 

where is the classical solution of □ 

Corollary 6.14. Assume (A^), g S e > 0, and Vg = B^[g]. 

(a) If {g)~ G A^(K^), then e /jj„(ue)" dx < /^^(s)" dx. 

(b) If g € Li(K^), then e ke| dx < /jj„ \g\ dx 

Proof, (a) Note that {g)~ e L^(IR^) implies that (— 5 )'*' € L^(R^). Since B^[—g\ = 
—B^[g\, we have by Proposition 6.13 that 

dx = £ f {B^[-g])+dx< [ {-g)+= f {g)~ dx. 

Jr" Jr" Jr" Jr" 


(b) Follows by noting that (xe)'*" + (x^) = |xe| 

Below, we collect the main results for (|3.1|). 


□ 


Theorem 6.15. Assume (A^ ), g & and Vg G L[q^(K'^) is a distributional 

solution of ( |3.1[ ). 

(a) If (g)^ G then 


: f (x£)+ dx< [ (g)+ dx. 

Jr" Jr" 


(b) If g > 0 a.e. on , then x^ > 0 a.e. on 


Proof, (a) Let cos G C° 
By assumption. 


be defined in (1.7), and let Vg^s = Vg-^ujs G 
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for all e C° 
equation 


. Taking ip{y) = ujs(x — y) for x € , we get the pointwise 


eVe,S - ^'^[Ve,s] = 9 * UJs 

Note that 0 < (g * < {g)~^ * cos & 

Proposition |6.13| gives 




(see e.g. Lemma 5.1 in [3S] ), so 


' dx < / (5 * ws)'*' dx. 

7b™ 7b™ 


Then by Patou’s lemma 


e / liminf(uE a)’*’dx < liminf e / (up 5)“''da; < liminf / ( 0 )“*" * dx. 

Since (•)+ is continuous, (g)+ € and the properties of 

mollifiers yields 


R™ 


(up)^ dx< (g)^ dx. 

7b™ 


(b) Note that —Vg solves (3.1) with right-hand side —g. If —g < 0 a.e. on , then 
(—g)+ = 0 e L^{R^). By part (a), we deduce that e/ rn(— da; < 0, and hence 


that —Up < 0 a.e. on ^ 


□ 


We are now ready to prove our main theorem for the elliptic equation (3.1). 


Proof of Theorem \3.1\ (a) By the assumptions and Proposition 6.8 (d), for every 
r > 0 , there exists a unique classical solution Vg^r € 

for all 


e||77“up,,|U~ < WD^gh^ 


of ( 6 . 1 ) satisfying 
a G N^. 


An Arzela-Ascoli argument as in the proof of Proposition 6.12 (in this case 
combined with a diagonal extraction argument), shows the existence of classical 


solutions Vg G C{)°( 


pN\ 


of (3.1) satisfying 


e|!77“up|Uoo < ||D“ 5 ||io= for all a G 
Moreover, Corollary |6.11| ensures that the classical solutions Vg are unique. 

(b) Existence of -solutions: Let S > 0, gs = g * uis G where ujs is 


defined by (1.7), and Vg^s G C'{)°(K'^) be the solution of (3.1) with gs as right hand 
side. By Remark |6.1| (a), a difference of solutions is also a solution, and then by 
Corollary [6T4l (b), 

£:\\vg,Si - Vg^S 2 \\L^ < \\ 9 S 1 - gs 2 \\L^ for every ( 5 i,( 52 > 0 . 

Hence, {up_ 5 } 5 >o is Cauchy and there exists Vg G L^(K^) such that ||up_ 5 —UpH^i —0 
as (5 —)■ 0 +. 

Since Vg^s satisfies (3.1) with right-hand side gs-, 

e I Vg^stpdx— I Vg^s^^[fp]dx = f gstpdx for all '0 G (^“(K^), 

7b™ ’ 7b™ ’ 7b™ 

and since Vg^s,gs —>■ Ve ,9 in as <5 —)■ 0 +, we send J —)■ 0 + and find that Vg is 

an L^-distributional solution of (3.1). 

Uniqueness: Note that C Consider two distributional solutions Vg,Vg 

of (|3.1|) with right-hand sides g,g G L^{R^). If g — g = 0 a.e., then Vg — Vg = 


Bg[g — g] =0 by Theorem 6.15 (b). 

-estimate: By the assu mpti ons, we can take Vg G L^(K^) C Aoc(®^) and 


g G L^(R^). Then Theorem 6.15 (a) gives 


e||(^s)+IUi<||(5)+||Li. 
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A similar argument as in the proof of Corollary |6. 14 concludes the proof. 

(c) Existence of L°°-solutions: Proposition | 6 . 8 | (b) ensures that there exists a unique 
a.e. solution Vg^r G of 

and e||we,r-||L“ < IIsIIl”- Then, by Alaoglu’s theorem there exists vf S 
such that, up to a subsequence, ^ wj in as r„ —)• 0+. That is, 

lim / Ue_r„'0da;= / vfipdx for all ip G 

*-7.^0+ Jrjv ’ Jrn 

To finish the existence proof, we need to show that vf is in fact a distributional so¬ 
lution of dnj). Consider a function 7 e then G 

(see Lemma 3.5 (b)). Since Ue,r„ is a pointwise a.e. solution and 1 ^ 5 ,t-„, \ve.r b G 

we have by integration and self-adjointness of (cf. Lemma 3.5 and 
Remark (b)) that 


£ [ i^e,r-„ 7 da;— f [ 7 ] dx = f g^dx for all 7 eC'“(M^). 

The weak* L°°-convergence ensures that 

/dx for all 7 GC“(M^). 


lim 

r-ji—>-0+ 


/ Xg_r„7dx= / Veji 

Jm." 


By Remark 6.1 (b), we have that, for any 7 G [ 7 ] —)• C^b\ L^{R^) 

as Tn -G 0+. Then, since \\ve,rJ\L'^ < jllslU”: we get as r„ -t 0+ 

Jr" 

= [ Xe,r„>C'"[ 7 ] dx-f [ - C^)[y]dx ^ [ v;£^b]dx, 

Jr" Jr" Jr" 

for all 7 G This shows that 

e [ zqqdx— f ZG£'^[ 7 ]dx= f 37 dx for all ^ G 

Jr" Jr" Jr" 

that is, vf is an L°°-distributional solution of ([ 33 . 

Uniqueness: Note that L°° C Lioc- Consider two distributional solutions Ve,Vg 
of (3.11 with right-hand sides g,g G li g — g = 0 on , then Vg — = 0 


by Theorem 6.15 (b). 

L°°-estimate: Observe that ±d||g||ioo G C are distributional 

solutions of (3.1) with ±||g ||L°° as right-hand sides. Moreover, —|| 5 ||l~ < g < 


ioo. Then Theorem 6.15 (b) gives |xe| < r|| 5 ||L=o. 


□ 


This section is concluded by a proof of the self-adjointness of Bp-. 


Proof of Lemma \3.f\ Let fs = f*0Js and gs = g * ujs where ujs is defined by ( 1.7|). 
Then fs G C^(W) C ^nd g G (M^), and then by Theorem 

(aHc), Bplfs] G n lT2.i(RAr)^ ^ and 


3.1 


sBPlfs] - Cf^[Bp[fs]] = fsix) in 
[ 55 ] - [Bp [ 55 ]] =gs{x) in 


oW 


By the regularity and integrability of the terms of the equations (cf. Lemma 3.5), 
we may multiply the first equation by Bp[gs] and the second by Bp[fs], and then 
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integrate both equations in x over By self-adjointness of (Lemma [3.5| (c)), 
we then find that 


/r« 


fsB^[gs] dx = / [fs] - [B^[fs]])B^ [g^] dx 


(6.7) 


/B« 


{eB^[gs]-C^mgs]])B^[fs]dx 


= [ 9sB^[fs]dx 

To pass to the limit as 5 —)■ 0^, we first subtract equations to find that 
eB^ [f] - eB^ [/,] - [B^^ [/] - B^ [fs]] = f - fs in V'{R^), 

and hence by Theorem 3.1 (b), linearity, and properties of mollifiers, 

e\\B^[f]-B^[fs]h^=e]]B^[f-fs]\\L^<]]f-fs\\L^^O as <5^0+. 
On the other hand, by Theorem 3.1 (b) and (c), and properties of the mollifiers, 

e\mfs]\\L^ < II/IIl^, e\mfs]\\L^ < ||/|k~, e\mgs]]\L^ < ||g|U~, 

and gs ^ g a-e. Using -convergence for the /-terms and the dominated conver¬ 
gence theorem for the g-terms, we may send 5 —>■ 0"^ in (6.7) to get the result. □ 

Appendix A. Technical results 
A. l. Proof of Liouville type of theorem. 

Proof of Theorem \3.SI[ By the definition of distributional solutions, 

[ v{y)C^^[f;]{y) dy = 0 for all G 

Jr" 

Let X G take i/'(y) = uJsix — y), whe re uj s is defined in (|1.7[ ), and let vs 
V *ujs G Co(M^) n C^(R^). By Lemma ^ , 

Fubini’s theorem to see that 


(b), C^lf)] G and we may use 


(A.l) 


£^[ 115 ] (x) = 0 for every x G 




Assume that there exists an x G such that vs(x) 7 ^ 0. We only consider the 
case vs(x) > 0; the proof in the other case is similar. Then M := sup^jg^iv vs > 0, 
and since vs G C'o(R^) there exists an xq such that 

0 < M = max vs = vs(xo). 


By equation (A.l) and Lemma 3.5 (b), we then find that 


0 = £^[ x 5 ]( xo ) = / (vs{xo + z) - vs{xo) - z ■ Dvsixo)] dg{z) 
J\Z\<K ^ ' 

+ (x5(xo-I-z) - X5(xo)) d/r(z) 

J\Z\>K ^ ' 




'■df,{z 


' \z\>K 


(^z; 5 (xo + z) - dg{z). 


Take any zq G supp y. By definition, zq ^ 0 and y,{B{zo,r)) > 0 for all r > 0. 
Hence we can take r, k G (0,1) small enough such that 

B{zo, r) n {0 G : |z| < k} = 0. 
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Since k < 1, vs{xo + z) — M <0, and B{zo,r) C {z G : \z\ > k}, the above 
inequality yields that 

[ (Mxo +z)-M)dfi{z)> -WD^vsW^^rg, [ |zpl|^|<«,d/4(z). 

JB{zo,r) ^ ^ ^ D J|^|<i 

Taking the limit as k —> O'*" using Lebesgue’s dominated convergence theorem (the 
integrand is dominated by jzp which is integrable by <B) gives 




(vs{xo + z) d/r(z) - Mn{B{zo, r))J dfi(z) > 0. 


Then by continuity, vs(xo + z) = vs(xo + zg) + A (|2 — zg\) in B{zo,r) for some 
modulus of continuity A, and we find that 

^^(cco + zq) + A(r) > ^ / vs{xo + z) dfj,{z) > M. 

H[B{zo,r)) Jb(zo,t) 

Hence, we may send r —0+ and get that vg{xo + zq) > M. It follows that 
vs{xo + Zq) = M since M is the maximum of vg. 

Repeating the above argument, we find that V 5 {xq + nzg) = M for every n S N, 
and thus 

limsupr! 5 (xo + nzg) > M > 0. 

n—>-oo 

This is a contradiction since lim|,j|_^oo 115 ( 0 ;) = 0. So, we conclude that vs(x) = 0 
for every x G 

By the properties of mollifiers, vg —?■ v locally uniformly in as <5 —O’*', and 
hence it follows that also v(x) = 0 for every x G □ 


A.2. Proof of a measure theory result. 

Proof of Lemma \6.J\ Remember that we defined 


Mi{A,B)= [ ( [ dii(a;)) dz = f v{B — z)dz, 

M2{A,B)= f ( [ dii(x)) dz = f i^{A — z)dz, 

JB \Ja-z j jB 

and that we want to show that Mi (A, B) = M2 (A, B). 

Consider the set C C R^^ defined as 

C = {{x,z)gR^^ :zGA, xGB-z}. 

Furthermore, define the sets 

S = {x = xb - xa'- XA & A, xb & B} = [J (B - xa), 

= {z G A : X G B — z} = {z G A : z G B — x} = An {B — x). 

Note that C can also be expressed as 

C = {(a;, z) e R^^ : a; e 5 , z € G^}. 

Then 

(A.2) 

Mi{A,B)=f (f lA{z)lB-z{x)di^{x)\ dz = f (f lc(x, z) dz^(a;)^ dz 
VJrn / JRN / 

= /" ( f '^c{x,z)dz\ dB{x) = [ ( [ ls(a;) lGjz)dz) dv{x) 
Ar« v./r« / ./r« \Jh« / 

= [ ( [ dz^ diz{x) = [ |G2;|di^(a:), 

Js \Jg^ j Js 
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Figure 1. 


where the third equality follows by Tonelli’s theorem (the tensor measure is a non¬ 
negative Radon measure), and | • | denotes the Lebesgue measure on 

We can proceed in the same way to change the order of integration in the ex¬ 
pression for M 2 (A, -B), but first we make use of the symmetry of v 

M 2 {A,B)= f i'{A — z)dz= f i'{—A + z)dz= [ ( [ du(a;)) dz. 

Jb JB Jb \Jz-A J 

Using the same technique we consider the sets, 

C = {(x,z) e : 2 e B, x£-A + z}, 

S = {x = xb + xa '■ XA & -A, xb & B} = 1^ (xb - A), 

xb^B 


Gx — ^ B \ X ^ —-A -\- ^ B z ^ A-\- x} — B C] (A -1- x). 


Second, we follow (A.21 to get 
(A.3) 


M2{A,B)= / |G,|du(x). 


Now, note that S = S. Moreover, Gx = A f] {B — x) is just a translation of 
Gx = B n (A -I- x). Then \Gx\ = \Gx\, since the Lebesgue measure is invariant 
under translations. (Consult Figure for a visual overview of the sets.) 

Finally, we can conclude by (|A.2[) and (A.3) that 


Mi(A,B)= / |G,|du(x)= / |G,|du(x) = M 2 (A,B), 

Js Js 

which completes the proof. 


□ 





































44 


F. DEL TESO, J. ENDAL, AND E. R. JAKOBSEN 


Acknowledgements. E. R. Jakobsen was supported by the grant Waves and Non¬ 
linear Phenomena (WaNP) from the Research Council of Norway. F. del Teso was 
supported by the FPU grant AP2010-1843 and the grants MTM2011-24696 and 
MTM2014-52240-P from the Spanish Ministry of Education, Culture and Sports, 
the BERC 2014-2017 program from the Basque Government, and BCAM Severe 
Ochoa excellence accreditation SEV-2013-0323 from Spanish Ministry of Economy 
and Competitiveness (MINECO). He would like to thank the Department of Math¬ 
ematical Sciences at NTNU for hosting him from September to December 2014. 
We also thank Juan Luis Vazquez for encouraging and supporting the stay of the 
second author at NTNU where this paper was started, and for useful comments 
and discussions. Finally we thank the anonymous referee for his suggestions and 
remarks that helped us improve our paper. 


References 

[1] N. Alibaud. Entropy formulation for fractal conservation laws. J. Evol. Equ,, 7(1):145-175, 
2007. 

[2] N. Alibaud and B. Andreianov. Non-uniqueness of weak solutions for the fractal Burgers 
equation. Ann. Inst, H. Poincare Anal, Non Lineaire, 27(4):997-1016, 2010. 

[3] N. Alibaud, S. Cifani and E. R. Jakobsen. Continuous dependence estimates for nonlinear 
fractional convection-diffusion equations. SIAM J. Math. Anal,, 44(2):603-632, 2012. 

[4] N. Alibaud, J. Endal and E. R. Jakobsen. On the duality between Hamilton-Jacobi- 
Bellman and degenerate parabolic equations. Work in progress. 

[5] B. Andreianov and M. Maliki. A note on uniqueness of entropy solutions to degenerate 
parabolic equations in NoDEA Nonlinear Differential Equations AppL, 17;109-118, 2010. 

[6] F. Andreu-Vaillo, J. M. Mazon, J. D. Rossi and J. J. Toledo-Melero. Nonlocal Diffu¬ 
sion Problems. Mathematical Surveys and Monographs, 165, AMS, Rhode Island, 2010. 

[7] D. Applebaum. Levy processes and Stochastic Calculus, Cambridge University Press, Cam¬ 
bridge, 2009. 

[8] B. Barrios, I. Peral, F. Soria and E. Valdinoci. A Widder’s type theorem for the heat 
equation with nonlocal diffusion. Arch. Ration. Mech. Anal., 213(2):629—650, 2014 

[9] J. Bertoin Levy processes. Cambridge Tracts in Mathematics, 121, Cambridge University 
Press, Cambridge, 1996. 

[10] P. Biler, C. Imbert and G. Karch. The nonlocal porous medium equation; Barenblatt 
profiles and other weak solutions. Arch. Ration. Mech. Anal., 215:497—529, 2015. 

[11] P. Biler, G. Karch and R. Monneau. Nonlinear diffusion of dislocation density and self¬ 
similar solutions. Comm. Math, Phys., 294(1):145-168, 2010. 

[12] M. Bonforte, a. Segatti and j. L. Vazquez. Non-existence and instantaneous extinction of 
solutions for singular nonlinear fractional diffusion equations. Calc. Var. Partial Differential 
Equations, 55(3):55—68, 2016. 

[13] M. Bonforte, Y. Sire and J. L. Vazquez. Existence, Uniqueness and Asymptotic behaviour 
for fractional porous medium equations on bounded domains. Discrete Contin. Dyn. Syst., 
35(12):5725-5767, 2015. 

[14] M. Bonforte, Y. Sire and J. L. Vazquez. Optimal Existence and Uniqueness Theory for 
the Fractional Heat Equation. Preprint, arXiv:1606.00873v2 [math.AP], 2016. 

[15] M. Bonforte and J. L. Vazquez. Quantitative Local and Global A Priori Estimates for 
Fractional Nonlinear Diffusion Equations. Adv. Math., 250:242-284, 2014. 

[16] M. Bonforte and J. L. Vazquez. A Priori Estimates for Fractional Nonlinear Degenerate 
Diffusion Equations on bounded domains. Arch. Ration. Mech. Anal., 218(l):317-362, 2015. 

[17] M. Bonforte and J. L. Vazquez. Fractional Nonlinear Degenerate Diffusion Equations on 
Bounded Domains Part I. Existence, Uniqueness and Upper Bounds. Nonlin. Anal. TMA., 
131:363-398, 2016. 

[18] H. Brezis and M. G. Crandall. Uniqueness of solutions of the initial-value problem for 
ut — /\ip{u) = 0. J. Math. Pures AppL, 58(2):153-163, 1979. 

[19] L. Caffarelli and J. L. Vazquez. Nonlinear porous medium flow with fractional potential 
pressure. Arch. Ration. Mech. Anal., 202(2):537—565, 2011. 

[20] J. Carrillo. Entropy Solutions for nonlinear degenerate problems. Arch. Ra¬ 
tion. Mech. Anal, 147(4):269-361, 1999. 

[21] E. Chasseigne and E. R. Jakobsen. On nonlocal quasilinear equations and their local limits. 
Preprint, arXiv:1503.06939vl [math.AP], 2015. 



NONLOCAL POROUS MEDIUM EQUATIONS 


45 


[22] S. CiFANi AND E. R. Jakobsen. Entropy formulation for degenerate fractional order 
convection-diffusion equations. Ann. Inst. H. Poincare Anal. Non Lineaire, 28(3):413-441, 
2011. 

[23] S. CiFANi AND E. R. Jakobsen. On numerical methods and error estimates for degenerate 
fractional convection-diffusion equations. Numer. Math., 127(3):447—483, 2014. 

[24] J. Droniou and C. Imbert. Fractal first order partial differential equations. Arch. Ra¬ 
tion. Mech. Anal., 182(2);299-331, 2006. 

[25] J. Endal and E. R. Jakobsen. Contraction for Bounded (Nonintegrable) Solutions of 
Degenerate Parabolic Equations. SIAM J. Math. Anal., 46(6):3957-3982, 2014. 

[26] L. C. Evans. Partial Differential Equations. Graduate studies in Mathematics, 19, AMS, 
Rhode Island, 2010. 

[27] G. Grillo, M. Muratori and F. Punzo. Fractional porous media equations: existence and 
uniqueness of weak solutions with measure data. Calc. Var. Partial Differential Equations, 
54(3):3303-3335, 2015. 

[28] M. A. Herrero and M. Pierre. The Gauchy problem for ut = when 0 < m < 1. 

Trans. Am. Math. Soc., 291(1):145—158, 1985. 

[29] H. Holden and N. H. Risebro. Front Tracking for Hyperbolic Conservation Laws. Applied 
Mathematical Sciences, 152, Springer, Berlin, 2007. 

[30] S. Kaliszewski, M. B. Landstad and J. QuiGG. Properness conditions for actions and 
coactions. Preprint, arXiv:1504.03394vl [math.OA], 2015. 

[31] S. N. Kruzkov. First order quasilinear equations with several independent variables. Math. 
Sb. (N.S.) (in Russian), 81(123):228-255, 1970. 

[32] O. A. Oleinik, A. S. Kalashnikov and Y.-l. Gzou. The Gauchy problem and boundary 
problems for equations of the type of non-stationary filtration. Izv. Akad. Nauk SSSR. Ser. 
Mat. (in Russian), 22:667-704, 1958. 

[33] E. Di Nezza, G. Palatucci and E. Valdinoci. Hitchhiker’s guide to the fractional Sobolev 
spaces. Bull. Sci. Math., 136(5):521-573, 2012. 

[34] R. Metzler and J. Klafter. The restaurant at the end of the random walk: recent de¬ 
velopments in the description of anomalous transport by fractional dynamics. J. Phys. A, 
37(31):R161-R208, 2004. 

[35] A. DE Pablo, F. Quiros and A. RodrIguez. Nonlocal filtration equation with rough kernels. 
Nonlin. Anal. TMA., 137, 402-425, 2016. 

[36] A. DE Pablo, F. Quiros, A. Rodriguez and J. L. Vazquez. A fractional porous medium 
equation. Adv. Math., 226(2):1378—1409, 2011. 

[37] A. DE Pablo, F. Quiros, A. Rodriguez and J. L. Vazquez. A general fractional porous 
medium equation. Comm. Pure Appl. Math., 65(9):1242-1284, 2012. 

[38] A. DE Pablo, F. Quiros, A. Rodriguez and J. L. Vazquez. Glassical solutions for a loga¬ 
rithmic fractional diffusion equation. J. Math. Pures Appl., 101(6):901-924, 2014. 

[39] W. SCHOUTENS. Levy Processes in Finance: Pricing Financial Derivatives. Wiley series in 
probability and statistics, Wiley, Ghichester, 2003. 

[40] D. Stan, F. del Teso and J. L. Vazquez. Finite and infinite speed of propagation for porous 
medium equations with fractional pressure. C. R. Math. Acad. Sci. Paris, 119:62-73, 2015. 

[41] D. Stan, F. del Teso and J. L. Vazquez. Finite and infinite speed of propagation for porous 
medium equations with nonlocal pressure. J. Differential Equations., 260(2):1154-1199, 2016. 

[42] F. DEL Teso. Finite difference method for a fractional porous medium equation. Calcolo, 
51:615-638, 2014. 

[43] F. DEL Teso and J. L. Vazquez. Finite difference method for a general fractional porous 
medium equation. Preprint, arXiv:1307.2474vl [math.AP], 2013. 

[44] A. Tychonov. Theoremes d’unicite pour I’equation de la chaleur. Mat. Sb., 42(2):199-216, 
1935 . 

[45] J. L. Vazquez. New self-similar solutions of the porous medium equation and the theory of 
solutions of changing sign. Nonlinear Anal., 15(10):931-942, 1990. 

[46] J. L. Vazquez. Nonlinear Diffusion with Fractional Laplacian Operators. In: Nonlinear 
Partial Differential Equations, Abel Symposia, 7:271-298, Springer-Verlag, Berlin, 2012. 

[47] J. L. Vazquez. Recent progress in the theory of Nonlinear Diffusion with Fractional Laplacian 
Operators. Discrete Contin. Dyn. Syst. Ser. S., 7(4), 857—885, 2014. 

[48] J. L. Vazquez. The porous medium equation. Mathematical theory. Oxford Math. Monogr., 
The Clarendon Press, Oxford University Press, Oxford, 2007. 

[49] W. WOYCZYNSKI. Levy processes in the physical sciences. In: Levy processes, 241-266, 
Birkhauser Boston, Boston, 2001. 



46 


F. DEL TESO, J. ENDAL, AND E. R. JAKOBSEN 


(F. del Teso) BASQUE Center for Applied Mathematics (BCAM), Bilbao, Spain 

E-mail address : felix.deltesoOuam.es 

URL : http://www.bcamath.org/es/people/fdelteso 

(J. Endal) Department of Mathematical Sciences, Norwegian University of Science 
AND Technology (NTNU), N-7491 Trondheim, Norway 
E-mail address : jorgen. endalSmath.ntnu.no 
URL : http: //www. math. ntnu. no/ ■ j orgeen/ 

(E. R. Jakobsen) Department of Mathematigal Sciences, Norwegian University of Sgi- 
ENCE and Teghnology (NTNU), N-7491 Trondheim, Norway 
E-mail address : erjamath.ntnu.no 
URL : http: //www. math. ntnu. no/'er j / 



